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XLV. The Effect of Discontinuous Space Charge on Rectifier Characteristics 


By G. G. WILEs 


University of the Witwatersrand, Johannesburg* 


[Received in revised form November 24, 1954] 


SUMMARY 


The paper is theoretical. An expression is derived for the local current 
density in the blocking direction near any impurity centre that happens 
to be situated at an optimum distance from the metal. Numerical ex- 
amples are given. The effects of these large local currents on the average 
current through the interface is evaluated in the case of a uniform dis- 
tribution of impurity centres. When the semi-conductor has a dielectric 
constant of 10 the error at room temperature caused by assuming the 
usual continuous distribution of space charge is found to be very small 
with 1016 impurity centres per cm?, but becomes appreciable at 107°. 
The case. of a random distribution is treated by a more approximate 
method but this is considered to be a more satisfactory method than that 
used by Schottky. The results are much the same as for the case of a 
uniform distribution. 


t 
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$1. Tae Continuous Spack-CHaRGH APPROXIMATION 


WHEN a suitable metal and semi-conductor are in contact the well-known 
Schottky potential barrier (Schottky 1939) is formed by the charged 
impurity centres in the semi-conductor. The space charge is often 
regarded as a continuous uniform distribution of charge. When this 
approximation is made we have a field whose intensity / near the metal 
is given by 
| F2=8nNe(V,,+V Jl real) 
where N is the number of impurity centres per cm, ¢ is the electronic 
charge, « is the dielectric constant of the semi-conductor, Va is the 
potential applied externally in the blocking direction, and eV, 18 approxi- 
mately the height of the barrier when no current is flowing. Throughout 
this paper the symbol Ff’ will denote the field defined by eqn. (1). Tf 
zero potential is chosen at the energy level of the bottom of the conduction 
band of the semi-conductor when this is isolated and uncharged, then the 
potential V at a distance x not too far from the metal is 
V=f3 eee er ( 2) 


as in fig. 1. The barrier reaches its maximum height at the interface 


20 and the potential there is V=0. 


* Communicated by the Author. 
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On the diode theory of Bethe (1942) the electrons are considered to 
have a sufficiently long mean free path to pass right across the barrier 
without suffering collisions. When eV,,>7, the current J per cm? is 
given by 

2J /ngevp=exp (—eV ,/kT)—1 ee in ee 


where & is Boltzmann’s constant and 7' is the absolute temperature. 
The symbol v, is the average velocity of an electron in one direction in a 
Maxwellian distribution and so is equal to (2k7'/7m)"/? where m is the mass 
of an electron. The symbol n, denotes the number of carriers per cm? in 
the conduction band at x=0. 


FERMI |SURFACE 


ME TAL|SEMI-CONDUCTOR 


The top of the barrier due to a continuous uniform distribution of space charge 
is at the potential V—0. The image face lowers the barrier to ZH). An 
impurity centre in either of the positions x—X lowers the barrier to ZF. 
An impurity centre at the optimum position X, lowers the barrier to 
H,. Two impurity centres near X, lower the barrier to F,. 


On the diffusion theory of Mott (1939) the electrons are considered to 
have a short mean free path so that they suffer many collisions while 
passing across the barrier. In this case the current density J is given by 


J/ngev,F =exp (—eV ,/kT)—1 oe aay tee 


where v, is the mobility of the electrons. 


§2. THe Errecr or LowERING THE PoTENTIAL BARRIER 


The image force produces in effect a lowering of the top of the barrier 
from V=0 to V=£, as in figs. 1 and 2. An impurity centre situated near 
the interface produces a further lowering in its vicinity from V=E£, to 
Veh. 

This paper is concerned with single carrier theories of rectification. 
On the diode theory the effect of lowering the top of the barrier from 0 to 
Ey is to increase the current density from the value J given in eqn. (3) 
to the value jy given by jo/J=exp (eH,/kT). The effect of lowering the 
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top of the barrier from H, to E is to increase the local current density 
from j, to 7 given by 

jljo=exp [eC(H—E, kT]. . . . . . . (8) 


On the diffusion theory the effect of lowering the top of the barrier 
from 0 to Ey is to increase the current density from the value J given in 
eqn. (4) to the value j, given by 


jo/J =kT[eF [ exp (—eV/kT) dx]-} ee ee 


where the integration extends across the whole barrier. Similarly the 
effect of lowering the top of the barrier from H, to E is to increase the 
local current density from j, to 7 given by 


j/J=kT[eF { exp (—eV/kT) dx]}-}. Saris OAT) 
Equation (6) is derived by dividing the more familiar equation 
jo exp (—eV/kT) dx=nyv,kT [exp (—eV,/kT)—1] 


by eqn. (4). Equation (7) is derived similarly. 


The potential barrier at points in a plane z—=constant. An impurity centre 
lies in this plane at the optimum distance X, from the metal interface 
r=): 


The integral in eqn. (6) across the barrier lowered by image force has 
been considered by Landsberg (1951) and by Billig (1951). However, the 
barrier referred to in eqn. (7) is of a more general shape and will sometimes 
have two maxima as seen in figs. | and 2. We can find the order of mag- 
nitude of the integral across this barrier by using V=H in the range 
esto 7b. and V=Fe from eqn. (2) in the range x=E/F to U= 0. 
Using this method of approximation in both (6) and (7) we can divide 
(7) by (6) and find 

j[jo=% ©XP [e(H—E,)/kT'| fees «(8) 
where a is the comparatively unimportant factor 


a=(eH ytkl (eH kT). 


) fay Pa adare atk maonitiuiqde 
We conclude from eqns. (5) and (8) that as far as orders of magnitude 
are concerned the effect of lowering a barrier from the general level 
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V=E, to a local level V=E£ in the vicinity of an impurity centre is to 
increase the current density from the general value j, to a local value j 
given by eqn. (5). 


§3. THE Errect oF ONE IMPURITY CENTRE 


Let a point be chosen on the interface between the metal and the semi- 
conductor so that the normal from this point does not pass close to any 
impurity centre. Then fig. 1 shows how the maximum of the barrier - 
along this normal is depressed by the image force from =0 to V=E, 
giving rise to a current density jy. Let another normal be chosen so 
as to pass through an impurity centre at a distance X from the interface. 
Figure 1 shows how the maximum of the barrier along this normal 
is depressed into two maxima, one on each side of the potential well at 
x=X. We denote the potential at the higher of the two maxima by F 
and the local current density there by j. 

There is an optimum position x= X, for the impurity centre which makes 
the levels of the two maxima the same. We denote this level as =F. 
This is the lowest level to which an impurity centre can depress the barrier. 
For if the impurity centre lies further from the metal then the barrier 
is higher than #, on the metal side of the well. If the impurity centre 
lies closer to the metal then the barrier is higher than #, on the side 
away from the metal. It follows that when the impurity centre lies at 
the optimum distance X, then the local current density has the largest 
value j, that is possible near a single impurity centre. Figure 1 indicates 
this situation and fig. 2 shows more detail. The interface is the plane 
«=0 and the barrier is shown for points in the plane z=constant that 
passes through the impurity centre. On this representation the effect 
of the impurity centre shows up as a saddle-shaped pass through the 
barrier with a well at the top of the pass. We shall now calculate the 
local current density j, due to electrons that make their way through the 
pass. 

By analogy with eqn. (5) or (8) the current density j, is given in terms 
of the current j, by 


Jiljo=exp [e(H,—E,)/kT']. eye ra eat )) 


To evaluate this we need expressions for H, and 1 To find #y we add 
a term e/4e x to the V of eqn. (2) so as to allow for the image force. The 
barrier has its maximum where V has its minimum, and this occurs at 
v==(e/4eF)1/, This gives at the maximum of the barrier the potential 


p= (el fe)t/?. a aay a ot eda ee a 
To find #, we consider the potential along the normal that passes 
through the impurity centre at X,. We regard the field near x=0 as 


consisting of two parts. Firstly there is the part 2e/eX,” due to the 
impurity centre at X, and its image at —X 1. Secondly there is the part 
F, due to all the other impurity centres and their images. In practical 
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cases F, is not very different from F'. The potential along the normal 
can now be written approximately as 


V=(2e/eX 2+ F,)u+e/4er, «<X,y, 
V=eje(x—X,)+F yx, «>Xy. 


In the case x<X, the last term is the contribution due to the image 
force. This approximation for V is reasonable near the minimum of V 
although it breaks down when x approaches either 0 or X,. In the case 
«x>X, the terms due to the image force and due to the image of the 
impurity centre have been neglected. By equating the minimum .of V 
when «<X, to the minimum of V when «>X, we find the optimum 
position X, to be X,=c,(e/4<F,)!/2 where c, is a solution of the quartic 
€;'+ 8c,3+ 12¢c,7=32 which gives the value c,=1-18. Since (e/4eF)1/2 is 
the value of x where the barrier has its maximum when a continuous 
distribution of charge is assumed, we see that the optimum position X , for 
a single impurity centre is just beyond this maximum as shown in fig. 2. 
Using this value of X, we find the value of the potential at the maxima 
to be 


Ey=(2+4e,)(eFy/e)¥2. ashe sees ang) 


We have derived the current j, near a single impurity centre lying at the 
optimum distance from the metal. If this centre is part of a uniform 
distribution of centres then the closest neighbours will be at a distance of 
about N~-1/3 or more. Hence we draw a sphere OA of volume 1/N that 
touches the interface at O with its diameter OA passing through the 
impurity centre at X,. We assume that all the other impurity centres 
can be represented by a continuous uniform distribution of space charge 
outside the sphere. This space charge and its image produce the field F’,. 
Since the radius of the sphere is (47.N/3)-/3 the appropriate value for F, is 


FF 2(efe\an.N/3)2 0 2 = eee) 


Some numerical values calculated for a uniform distribution with 
-=<10 at 17°c are shown in table 1. The first column gives the concen- 
tration of the impurity centres, NV. The second gives approximately the 
applied voltage V,. The third shows the optimum distance X, where 
the single impurity centre is lying so as to produce its maximum effect. 
This column also gives an idea of the size of the pass since the length of 
the pass across the well from one maximum to the other is about 2-3X,. 
The width of the pass at the same energy level is much the same. The 
fourth column shows the depth of the pass from the crest of the barrier 
E, down to the mouth of the well #,. This is about 1-64). The last 
column shows j;/j). Here we see the tremendous increase 1n the current 
density that can occur in the immediate neighbourhood it an impurity 
centre when this happens to lie at its optimum distance from the metal. 

The values in table 1 refer of course only to the case €= 10 at 17°c. 
Values of j,/7) at other values of « and at other temperatures can, however, 
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be inferred from table 1. For the value of j,/j) remains practically the 
same under different conditions provided N(V,,+V,)/7* remains the 
same. For example the value j,/j,=6200 given for V,,+V,—50 volts 


will also be the value for V,,+V,=5 volts if N=10!8 or if «=4-6 or if 


T—163°K. 
Table I 
N Vota d©, Depth of pass ldo 

(per em?) (volts) (angstroms) | (electron volts) — 
1016 l 30 0-07 a | 
97 8) 0-24 340 . 

101 1 17; 0-13 22, 

20 8 0-28 990 

50 6 0-36 6200 

1018 1 10 0-22 200 

12 5 0-44 40000 

22 4 0-51 270000 


§4. THe Errect or A UNIFORM DISTRIBUTION OF IMPURITY CENTRES. 


In view of the very large current density j, that can be obtained near 
an impurity centre, the question arises as to whether the average current 
density over a large area of interface is anywhere near equal to the current 
density jg calculated on the usual assumption of a continuous uniform 
distribution of space charge. In this section we shall examine the average 
current density j,, when the impurity centres are distributed uniformly 
through the semi-conductor. In § 5 we shall examine the average current 
density j, when the impurity centres are distributed at random. 

With either of these distributions we see from eqn. (5) that the ratio j,./j5 
or j,/Jo 18 equal to 


average value of exp (eH/kT7’) ; 
(13) 
exp (€H/kT) 

where # is the minimum value of V along the normal drawn at each surface 
element of the interface. With either distribution we can choose any 
arbitrary point O on the interface and draw a sphere OA of volume 
1/N to touch the interface at O (see figs. 3 and 4). We can then assume 
that as far as potentials in the neighbourhood of O are concerned, all 
the impurity centres outside the sphere can be treated as a continuous 
distribution of charge outside the sphere. Inside the sphere there may be 
no impurity centre, or one, or more than one. If the point O is allowed 


to move to many different points on the interface, the average number of 


impurity centres in the sphere will be one. 
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In this section we shall assume that the continuous distribution of 
charge outside the sphere OA is uniform so that it produces near the 
interface a constant field F', given by eqn. (12). We shall find the average 
current density j, through those parts of the interface for which the 
number of impurity centres inside the sphere OA is one. We shall call 


ju the average current density given by a uniform distribution of impurity 
centres for the following reason. 


Fig. 3 


<—METAL INTERFACE 
The sphere OA contains on the average one impurity centre. If this impurity 


centre lies anywhere on the surface XB it produces a barrier of the 
same height along the line OA. 


Fig. 4 


<METAL INTERFACE 


OA is the same sphere as in fig. 3. If an impurity centre lies anywhere on 
the surface XC it produces a barrier of the same height along the line 
OA. A centre in the shaded region does not affect the barrier height 


along OA appreciably. 
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The first approximation in existing theories is based on the assumption 
of a continuous uniform distribution of space charge. A better approxi- 
mation is obtained by allowing the charge e distributed through the 
sphere OA to condense to a point charge e somewhere within the sphere. 
The limitations implied in this approximation are (a) that the field due 
to distant centres is constant over the interface, and (6) that only one 
centre contributes to the formation of each pass across the barrier. 
Now these are just the limitations that would be imposed by a distri- 
bution of centres that is uniform instead of random. For a distribution 
of points may be regarded as uniform when (a) there are no large scale 
fluctuations in the concentration of the points, and (6) no two points are 
separated by a distance of less than about N~V?. 

It is assumed in the calculation of j,, that the probability of occupation 
by the impurity centre is the same for all points within the sphere OA. 
This is equivalent to a third limitation on the distribution of the centres, 
for it excludes distributions in which the centres are arranged in a regular 
pattern. If the centres were arranged in a regular pattern there would 
generally be an unduly large number or else an unduly small number of 
centres lying in the critical region near the plane x= X,. 

In the calculation of j,,/j, from the expression (13) the average value of 
exp (eH/kT) is 

N f exp(eH/kT) dv See Sh ee Repke Ls 


where # is the minimum value of V along the diameter OA when the 
impurity centre lies within the volume element dv. The integration 
extends through the sphere OA and the factor V comes in as the reciprocal 
of the volume of the sphere. To perform the integration it is convenient 
to consider firstly the part of the sphere beyond the plane «=X, as in 
fig. 3, and secondly the part between the interface and the plane x= X, 
as in fig. 4. 

Firstly, then, the impurity centre lies in the sphere to the right of the 
plane «=X, as in fig. 3. When it lies on the line OA at a point =X, then 
the field near O is 2e/eX?+-F, so that the potential is 


V = (2e/eX?+ F,)x+e/4ex. 
The minimum value of this is found with the help of the expression above 
for X, to be 
H=(1+-84 57/0 2X9) e(ehie)\ 2 eee) 
As a check we can use the quartic defining c, to find that when X=X,, 
E reduces to EL, given by (11) as it should. 

When the impurity centre lies off the line OA at a distance 7 from OA 
it still gives with its image the same field at O provided it lies on the 
surface r?=X4/32?/3—a?, This surface is very similar to the sphere 
r= Xx—2?, both being surfaces of revolution about the line OA cutting 
OA normally at x«=0 and «=X. Hence when the impurity centre lies 
to the right of the plane x=X, between two neighbouring spheres of 
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diameters X and X+dX then E is given by eqn. (15). The volume 
element between the spheres is 


dvu=(n/2)(X2—X,?) dX. 
‘The integration goes from X=X, to X=OA. This contribution to the 
integral in (14) was evaluated numerically. 
Secondly the impurity centre lies in the sphere OA to the left of the 


plane «=X, in fig. 4. When it lies on the line OA at any point =X we 
assume for the maximum of the barrier the linear relation 

B=[1+(1+4e)X/XMeF ie? . (16) 
‘This relation was obtained from the following considerations. When 
X=X,, H must reduce to the value H, given by eqn. (11). When X¥—0, 
the impurity centre and its image will give no contribution to the maximum 
of the barrier. Since the field is then 7, # must reduce to (eF',/e)#/2 by 
analogy with eqn. (10). At intermediate values of X the assumption of 
a linear variation was found to be a reasonable approximation in particular 
numerical cases. 

When the impurity centre lies off the line OA we take the locus of its 
positions that maintain the value of # in (16) to be the surface of a sphere 
‘that cuts OA at «=X and that has its centre at D in fig. 4. The point 
D is at x=(1+2/c,)X, so as to coincide with the second maximum of 
the barrier when the impurity centre is on OA at ~=X,. When X is 
near X, this locus is reasonable because the maximum of the barrier then 
-occurs near D. Wherever the impurity centre may be on the locus, the 
‘potential at D remains at the value e/eR+F,(1+2/c,)X, where £& is the 
radius of the locus. When X is not near X,, any locus is reasonable since 
the impurity centre in such positions does not depress the barrier height 
along OA appreciably. The volume element between two neighbouring 
spheres that cut OA at X and X+dX is 

du=2n{X?—(2-+2e,)X,X+(1+2/c,)X 7] dX. 
The integration goes from X=0 to X=X,. This contribution to the 
‘integral in (14) was integrated exactly. 

To complete the integral in (14) we still have to consider the remaining 
volume that is shaded in fig. 4. When the impurity centre lies in this 
region its depression of the barrier along OA is negligible so that EH has 
the constant value (e/’,/e)"?. The volume of the shaded region is 

n(3/4aN )V3X ,2@—a( 1+ 2/¢,)X,°. 
This gives the final contribution to the integral in (14) for substitution 
‘into (18). 

Some numerical values for the case «-=10 at 17°c are shown in table 2. 
The values of j,,/j) are given in the third column, ),, being the current 
density:in the blocking direction averaged over a large area of the inter- 
face when the distribution of centres is uniform, and j, being the corres- 
ponding current density calculated on the usual assumption ofa eens 
uniform distribution of space change. The values close to 1! show that 
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when there are 1016 or 1017 impurity centres per cm® the approximation 
based on a continuous distribution of space change is a good one. But 
at a concentration of 1018 the approximation is not so good. At first 
sight one might expect the higher concentrations to behave more like a 
continuous distribution than the lower concentrations. But the revérse 
is seen to be true. The reason is that the higher concentrations provide 
more impurity centres in the critical region near the plane x=X, than 
the lower concentrations do. 

A few values of j,,/j) in table 2 are less than 1. This may be an error 
due to the approximations used. On the other hand it seems possible 
in principle for j, to be less than jy. For the effect of averaging 
exp (eH/kT’) as in eqn. (13) instead of averaging by the use of a continuous. 
space change as in § 1 tends-to make j,, exceed jy. But the effect of intro- 
ducing the image force before averaging as in eqn. (15) instead of intro- 
ducing it afterwards as in § 2 tends to make j,, less than jo. 

The values in table 2 refer only to the case «=10 at 17°c. In the case 
of j,/jq in table 1 we pointed out that the tabulated values can be adapted 
for other values of « and 7’. Unfortunately this does not seem possible in 
the case of j,,/j9. It seems that a separate numerical integration is required 
for every case. 


Table 2 
N Vos 2 Ju Jo Ju Jo dr Jo } 
(per cm?) (volts) accurate approximate approximate 

1016 1 0:99 1-00 1-00 

97 0:99 1-00 1-00 
1017 1 0:98 1-00 1-00 

20 1-00 1-00 1-01 : 

50 1:02 1-01 1-02 
1018 1 1-02 1-02 1-03 

12 1:63 1-49 1-67 

22 2-91 3:04 3:70 


§5. Toe Errecr or a Ranpom Disrripution or Impurtry CENTRES. 

This problem was examined in § 10 of a paper by Schottky (1942). 
He discussed first the effect of a series of impurity centres arranged in a 
line. They were sufficiently close together to allow electrons to pass from 
one to the other through the barrier. He concluded that this arrangement 
of impurity centres would occur so seldom in a random distribution that 
it would have no appreciable effect on the general current. 

He discussed next the effect of random fluctuations in the concentration 
of the impurity centres. A local increase in the concentration has the- 
effect of reducing the width of the barrier in this locality. If the con- 
centration increases everywhere by a factor 7? then the width of the- 
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barrier decreases everywhere by a factor 1 /r. Accordingly Schottky 
assumed that a local increase in the concentration in the ratio r2 decreases 
the local width of the barrier from the general width b to the local width 
_ bir. He took the local region to be that part of the barrier that is enclosed 
in a cylinder having its axis normal to the interface and having an area 
of cross section equal to 6?. At the next stage of the argument a smaller 
cylinder with an area of cross section (b/r2) is chosen within the original 
cylinder. A random increase in the concentration in the second cylinder 
in the ratio r? reduces the barrier width from the general value b/r over 
the first cylinder to the local value b/r? over the second cylinder. This 
argument is repeated for successive stages and at each stage a new 
cylinder is introduced within the previous one. After n stages the local 
width of the barrier has been reduced to b/r” over an area (Oir2= aay We 
arrive in effect at a cone with its base of area b? in the plane x= and its 
apex in the interface. The random fluctuations in the concentration of 
impurity centres within the cone reduces the local width of the barrier 
indefinitely until there is a break through of electrons at the apex due 
to the tunnel effect. The assumption has been made that the field near 
the apex of the cone is mainly due to the space charge within the cone. 
On the basis of a continuous uniform distribution of charge about 87% of 
this field is caused by charge outside the cone. Only the remaining 13% 
is affected by the increase of concentration within the cone. Therefore 
the decrease in the local width of the barrier has been overestimated. 

Schottky used as an example the case r=3. He found it necessary to 
carry the process through about twenty stages in order to get a part 
of the barrier thin enough for a break through by the tunnel effect. He 
showed how to calculate the probability of such a fluctuation in the con- 
centration on the basis of a Poisson distribution. In the earlier stages 
the cylinders are large enough to contain many impurity centres. But 
in the later stages the cylinders become too small for this condition to 
hold. For example in the seventeenth cylinder the width of the barrier 
is reduced from 35 to 28 angstroms. This is caused by a random 
fluctuation that increases the number of impurity centres in the cylinder 
from 0-36 to 0-56. It becomes clear that there is a need for an analysis 
that deals with whole impurity centres. 

In the earlier sections of this paper we used a sphere OA drawn at any 
point O on the interface, with a volume 1/NV where N is the number of 
impurity centres per cm*. All the centres lying outside this sphere can be 
represented by a continuous distribution of space charge. In the case of 
a uniform distribution of impurity centres we assumed that there would 
be in effect only one centre within the sphere. In the case of a random 
distribution of impurity centres the sphere will often contain no impurit Vv 
centre, often one, sometimes two, and occasionally three or more. When 
there are n impurity centres in the sphere they will be most effective when 
they are close together and acting like a single centre of charge ne at an 
optimum distance x=X,,. They will then depress the pass through the 
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barrier to a potential V=Z,,, and give rise to a local current density j ,- 
By generalizing the previous theory we find eqns. (9) and (11) replaced by 


J nlJo= XP [e(# n—Ho) [kT] 7 7 i (17). 
where E ,=(2n1/2+-4¢,,)(eF,/€)¥/. 5s, 4 


Here c, is a solution of the quartic c,4+8,,1/c,3+4(4n—l)c ,?=32n. 
giving the values shown in table 3. 


Table 3 
n il 2 3 4 > 50a 
is : 1:18 1-22 1-24 1-25 i 47 


Incidentally the optimum distance X,, has the value 
Apa teaer 
so that the optimum position for a multiple centre is practically the 
same as the optimum position for a single centre. 

The average current density j,, has been given in eqn. (13). The sphere 
OA can be used as before to find the average value of exp (eH#/kT’). There 
will be n impurity centres in the sphere with a relative frequency of 1/n!. 
When n=0, H=(eF,/e)1/? with F, given by (12). When n=1 we average 
by expression (14). When n=2 the integral (14) is replaced by a double 
integral, when n=3 by a triple integral, and so on. But this process 
involves the numerical integration of multiple integrals. In order tc avoid 
the long computations this accurate method was abandoned and the 
following approximate method was used instead. 

The comparatively large sphere OA was replaced by a small critical 
region near the optimum position X, opposite the fixed point O. A sui- 
table volume for this region was found by trial and error to be v= X ,3/10. 
The average current density at O is obtained by assuming that when n 
impurity centres lie within the critical region the current density is j,, 
given by (17) and (18), and when no impurity centres lie within the critical 
region the current density is simply jp. 

The method has been tested on the case of a uniform distribution of 
impurity centres with «=10 at 17°c. In this case there cannot be more 
than one impurity centre in the critical region and the probability of 
finding even one there is Nv. The appropriate value for F, is given by 
eqn. (12) and the values of j ,,/j) obtained by this method are shown in the 
fourth column of table 2. They compare well with the accurate values 
in the third column. 

When the method is applied to the case of a random distribution of 
impurity centres we have the number of centres in the critical region 
governed by a Poisson distribution. The average number of impurity 
centres in regions of volume v is Nv. Hence the probability of finding n 
impurity centres in the critical region is (N"v"/n!)exp(—Nv). The 
appropriate value for F, is F,==F since the volume v is too small for a 
correction of the type in eqn. (12) to be necessary. The values obtained 
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for j,/jo with e=10 at 17°c are given in the last column of table 2. It is 
seen that the ratios are a little higher than they are in the case of a 
uniform distribution, as we might expect. 

It was demonstrated above that v—X,3/10 is a suitable volume for the 
critical region when this contains only one impurity centre. The question 
arises as to whether this is a suitable volume when there is more than one 
impurity centre in the critical region. The answer seems to be that 
it does not matter whether it is the correct volume or not in these cases. 
For when j,/j) differs appreciably from j,,/j, an examination of the 
numerical work shows that the difference between the two is not due to 
clusters of impurity centres in the critical region. The clusters occur too 
seldom. Even double impurity centres account for only about 1° of 
the difference. The remaining 99° of the difference is due to the use 
of F,=F in one case and the use of F, corrected by eqn. (12) in the 
other case. Physically this means that when a single impurity centre 
lies in the critical region, which occurs comparatively often, its effect 
is enhanced when other centres are allowed to lie near it at random. 
In the pictorial representation of fig. 2 these other centres lower the 
pass indirectly. The pass is still caused by the single impurity centre 
in the critical region but the other centres lower the barrier in which the 
pass is made. 

No allowance has been made in this paper for electrons that break 
through the barrier by the tunnel effect. When the distance from the 
metal to an impurity centre is a little greater than X,, the barrier is 
higher on the metal side of the well than it is on the side away from the 
metal. In cases when X, is small (see table 1, column 3) the barrier on 
the metal side may be narrow enough for electrons to break through by 
the tunnel effect. Then the electrons contributing to the current need 
only sufficient energy to surmount the lower part of the barrier on the 
side away from the metal. This will give rise to local current densities 
greater than those indicated by j,/jo in table 1, and to average current 
densities greater than those indicated by j,,/j) and j,/jo in table 2. 
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XLVI. A Contribution to the Theory of Nuclei with a Diffuse Surface 


By Z. JANKovIC* 
CERN, + Theoretical Study Group { 


[Received November 23, 1954] 


Tue theory of neutron reactions has recently made significant advances 
by introduction of the complex square-well potential (Feshbach, Porter 
and Weisskopf 1953). In this connection it seems worth while to investi- 
gate the consequences of replacing the sharp nuclear surface by a potential 
which varies more gradually. The aim of this letter is to enumerate 
only some results which the author has obtained while a more systematic 
and detailed exposition will be given in a separate article. 

The nuclear potential is assumed to vary continuously over a range, 
5, from its exterior value V=0 to the interior value 


V=—V,(1+76€). rome a Ne: 
For simplicity a linear variation of V in the surface region has been 
assumed (cf. fig. 1). From the solution of the Schroedinger equation the 
logarithmic derivative f and its derivative with respect to the energy 
have been determined. As in the region of low energies the s-neutrons 
play the most important part in nuclear reactions, we quote the logarithmic 
derivative for this case : 


FF+GG)—é/3(#G—FG 
Regia Se 


es (2) 
(FG—FG@)+é/3(FF+G4G 
In(\=(R43) | 
with P=—QJ 1/3(@) +p EJ 1/3'(@)+ £9 4/3" (@) | Or 45 


G= pI 1/3(4) +S 13 (4) +489 13'(@), J 
where J 4,3 is the Bessel function of order one third and the barred quanti- 
ties contain the Bessel functions of order minus one third. 

cot y/(V)) Ri coth? b— 1] 1 
coth b[cot?\/(V,))R+1] é 
cot?4/(V>)R-+coth? bh: © 48./V,’ 
Shoe 2) 
(= BVM ey D=E2V (VIB, O=BYVo = 3S Tee gps 6) 


(€ is assumed to be a small quantity compared with unity). 
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On the basis of these expressions, I" nlD, Sse) Ty, G, have been calculated 
by means of the well-known procedures (Blatt and Weisskopf 1952, 
Feshbach, Porter and Weisskopf, published 1954). The results of these 
_ calculations are summarized in the following diagrams, which plot the 
various quantities as a function of the radius R; the calculated cases 
have been chosen in such a way as to study the effect of independent 
variations of the several parameters occurring in the theory. 


Fig. 1 


The assumed shape of the nuclear potential. 


The data suggest that four maxima (Adair 1954) of the various cross 
sections occur in the range of the periodic system of the elements, i.e. 
between 0<R<10, and that the fourth maximum is important for the 
description of the behaviour of the heavier nuclei. As could be seen 
from the theory, the parameter €(€<1) does not influence appreciably 
the position of the maxima in this range, and the Ryax Values (determined 
by the condition that for R=R,,x, Re(f) vanishes and —Im/(f) is near 
its minimum value) could be determined from the requirement 


1 J’, (a) 
cob/(Vo)Ruax=— 35 7V- ~ Traala)’ oad () 


These Riyax Values compared with the corresponding values for the case 
au . . 
5—0 show a shift to smaller R values, the difference being greater than 
$/2 and increasing with V, and 6. Ae 
fe comparing the corresponding diagrams the main qualitative effects 
of changing only one of the parameters V,, 6, € may be obtained : 
(a) V,: greater V, shifts the curves to the smaller values, the distance 
hereon maxima becoming smaller and the width is decreased ; the 
i inima i nall one. 
influence on the value of maxima and minima Is a SI | 
- 6) 8: greater 5 shifts the curves to the smaller # values : all the cross 
sections are increased as is also the widths of the ‘resonances aa 
(c) €: inerease of € does not shift the curves in the first approximation ; 


values near the minima are increased, but maxima are depressed and the 


width is enlarged. 
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Fig. 2 
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The calculated total cross section for incident s-wave neutrons is plotted as a 
function of the radius parameter R (in units of 10-13 em). The curves 
refer to a neutron energy which varies with R so that kR=0-2, where 
k is the wave number of the incident neutrons. The parameters of the 
potential are for (a) the dashed curve, V,—42 mev, 85=0, £=0-03 ; the 
solid curve, V,—42 mev, 8=2 x 10-13 cm, €=0-03 (the curve shifted for 
6/2 to the right); for (6) the dashed curve, V)=—30 Mev, 5=2 10-13 cm, 
€=0-05; the solid curve, V,=30 Mev, 6=2 x 10-8 om, €=0-03. 


ak ee 


(2) 


(6) 


The ratio of neutron width I’, (reduced to 1 ev neutron energy) to the level 
spacing D is plotted as a function of the radius parameter # (in units 
of 10-18 cm). The parameters of the potentials used are the same as 


for the corresponding curves in figs. 2 (a) and 2 (5). 
2D4 
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In the diagrams, R was intentionally left as the independent variable, 
because of the several possible relations between this parameter and the 
atomic number A. Perhaps the best definition of the nuclear radius p in 
our case would be 

p=R-+k ; (32 Ras ete  ee 
(x is a constant between zero and 6) and then assume p (or &) to be 
proportional to A1/3 
plor h)=r, Ay: Er Eee 
In either case we may remark that the maxima are at equal distances 
7//V, from each other (irrespective of 5), i.e. 
R,—B,_=16(A A8—Aj91)=2/V Vo. il a oe 
Thus, from A, and A, we can infer A, and also the value for 79./ Vo. 
Now, if either r, or V, is given, R; can be determined from (8) and then 6 
from (6). The remaining parameter, £, should be adjusted to yield the 
experimentally observed peak to valley ratios of the cross sections (Feshbach, 
Porter and Weisskopf, published 1954). 

Finally we may remark that the inclusion of the nuclear deformation 
may also be important in the analysis of the experimental data (A. Bohr 
and B. R. Mottelson 1953). 

The author wishes to express his appreciation of the arrangements which 
made possible his participation in the CERN Theoretical Study Group 
in Copenhagen. He is very grateful to Dr. B. R. Mottelson and to 
Dr. A. Bohr for raising this question and for their interest regarding this 
work. 
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XLVI. The Disintegration of 7Be by Slow Neutrons 


By R. C. Hanna 
Atomic Energy Research Establishment, Harwell t 


[Revised MS. received November 26, 1954] 


SUMMARY 

The reaction 7Be (n p) 7Li has been identified for neutrons incident with 
thermal energies. The cross section has been measured as 5:3-40-8 x 10-29 
em?. The cross section for the reaction 7Be (n «) 4He does not exceed 
1 barn. The results show the parity of 7Be to be odd. 

The activity of 7Be source has been diminished by neutron bombard- 
ment. ‘The cross section for all processes leading to the destruction of 
*Be is found to the 4:8-L0-9 x 10-29 em2. 

A comparison of the 7Be (n p)’Li and 7Li(p n)7Be cross sections gives 
evidence for a level in *Be near the neutron threshold of width not 
exceeding 30 kev. 


§ 1. INTRODUCTION 
*Be is a K capture body of half life 53 days; in 11% of disintegrations 
the daughter nucleus, ‘Li, is left in an excited state at 478 kev and 
y-radiation is emitted (Turner 1949, Richards and Williamson 1949). 
From the energy balances of the reactions 
*Li+p=’Be-+n—1-65 Mev 
*Li+p=*Be+y-+17-2 Mev 
8Be=2a+80 kev 
the following reaction may be expected to occur when “Be is bombarded 
with thermal neutrons : 
7Be-+-n=‘Li+p+1-65 Mev 
7Be+-n=7Li*+p+1.17 Mev 
7Be-+n=2«-+18-9 Mev. 

The mode of disintegration of the excited *Be compound nucleus is 
determined by its parity. A *%Be nucleus of even parity and spin can split 
into two « particles, and the energy release for this reaction is so much 
greater than for proton emission that it would be expected to predominate. 
A ®Be nucleus of odd parity will de-excite by proton emission ; competition 
from y-ray emission would be expected to be small. Cross sections for 
neutrons of thermal energies are appreciable only for the incident neutron 


+ Communicated by the Author. 
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wave of zero angular momentum. Thus the parity of the compound 
8Be nucleus will be the same as that of the ground state of 7Be. A de- 
termination of the relative magnitudes of the cross sections for the 
np and nz reactions thus leads to an assignment of parity to the ground 
state of “Be. 

The measurement of the proportion of protons which leave the “Li 
nucleus excited is not readily interpreted (cf. 1°B (n «) 7Li). 

From the reciprocity theorem the cross sections for the reactions 
7Be (np) 7Li and 7Li (pn) 7Be for the same excitation of the compound 
nucleus are in the ratio 


o(n p)/o(p n)=P,?(2S y+ 1)(2J,+- 1)/P,?(2J ge+ 1)(2J,+ 1). 


In principle, therefore, the spin of 7Be can be determined. From 
measurements of the cross sections for the 7Be (n«) and 7Be (np) reactions 
both the spin and parity of the 7Be nucleus in its ground state can be 
deduced. 


§ 2. SruprES on 7Be+-n BY PARTICLE COUNTING 
2.1. Preparation and Assaying of *7Be Samples 

7Be was prepared by bombardment of a block of pile graphite with the 
175 Mev internal proton beam of the Harwell cyclotron. After chemical 
separation the *Be was deposited on platinum sheet as a ‘ weightless’ 
source. 

The mass of *Be was determined from the counting rate due to 478 
kev y-rays detected in an organic scintillation crystal relative to the rate 
due to 411 kev y-rays from a source of Au; the source strength of the 
latter was determined by coincidence counting. Corrections were made 
for internal conversion and the presence of harder radiation in the spectrum 
of Au (Cavanagh 1951). The difference of crystal detection efficiency 
for Au and “Be y-rays was calculated from the Klein—Nishina formula. 

Pulse height spectra taken with a sodium iodide detector, and half-life 
measurements, showed the 7Be source to be radioactively pure. The 
absolute activity was found to be 5-4-+-0-1x10® y-rays/second, which 
corresponds to a mass of 3-76---0-40x10-® g 7Be. The error is largely 
due to the 10° uncertainty in the value of the branching ratio in 7Be 
decay. 


2.2. The *Be (np) Reaction 
2.2.1. Method 


The *Be was mounted on the source electrode of a parallel plate gridded 
ionization chamber, and bombarded with a collimated beam of thermal 
neutrons from the Harwell pile B.E.P.O. A linear amplifier and four 
channel pulse amplitude analyser were used to study the spectrum of 
pulses from electron collection in the chamber. 

With the very small mass of target material available, the observation 
of disintegration particles of the expected energy was necessary but not 
sufficient evidence for the occurrence of the 7Be (np) reaction. At 
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sufficiently reduced pressures of the chamber filling gas, the reaction pro- 
ducts pass out of the collecting volume of the chamber and the associated 
pulses are reduced in height. From a study of pulse height spectra as a 
function of gas filling pressure, both the energy and the range of the 
particles responsible for a peak in the spectrum could be found, and thus 
the type of particle identified. 

A definite identification of the 7Be (np) reaction was made from the 
following observations. 

I. The spectrum of pulses showed a peak at an energy of 1:44 Mev 
when the *Be was present in the chamber. In its absence there was no 
peak. 

2. The variation in the shape of the spectrum as the gas filling pressure 
was reduced showed that the peak was due to protons. 


Kigml 


Counts 
per 
channel 


500 


Volts 


Pulse height spectra from the reaction BM(n a)" La. 
o 76 em of mercury filling pressure 
x 11 cm of mercury filling pressure 


3. Cadmium difference curves showed that the peak was due to 
incident thermal neutrons. . 

4. Since the proton energy was greater than the energy corresponding 
to the neutron proton mass difference (0-78 Mev) the target nucleus must 
have been heavier than the product nucleusi.e. unstable against K-capture. 
7Be is the only neutron deficient product of the proton bombardment of 
carbon to give the correct energy balance. 
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5. The intensity of the proton group was measured over a period of 
time. After normalizing to constant incident neutron flux, the intensity 
of the proton group was found to decrease with the expected 53 day 
half-life. This confirms the instability of the target material ; the value 
of the half-life is further proof that the reaction studied is due to “Be. 


2.2.2. Results 


In fig. 1 are shown the pulse height distribution curves obtained when 
a thin sample of 1°B was bombarded with thermal neutrons, for gas 
filling pressures of 76 cm and 11 cm of mercury. The voltages applied 
to the chamber electrodes were adjusted to be proportional to pressure. 
A small shift occurs in pulse height; the main features are seen to be 
independent of pressure. 


Fig. 2 
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In fig. 2 are shown the pulse height distribution curves obtained when 
"Be was bombarded. Curve 1 was taken at 76 cm pressure when the 
range of protons of the peak energy is less than the minimum chamber 
dimension ; curve 2, taken at 51 cm shows a reduction in the peak 
counting rate. Very few particles lose their full energy in the collecting 
volume at 35 cm pressure (curve 3). At 11 cm pressure (curve 4) a con- 
tinuum remains, believed due to «-particles from the !°B (n «) reaction 
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emitted from depth in the chamber electrodes and walls. The collecting 
volume of the chamber is not well defined ; the shape of the pulse height 
spectrum from a group of monoenergetic particles of range greater than 
the minimum dimensions of the collecting volume cannot be calculated. 
It may be noted that in the geometry used the calculated value of the 
pressure required to stop a proton of energy 1-44 mev from reaching the 
grid was 58 cm of mercury. The corresponding pressures for deuterons 
tritons and «-particles of the same energy were calculated to be 37, 30 
and 10 cm respectively. At 35 cm some protons could lose their whole 
energy within the space bounded by the source and grid electrodes. 
At 11 cm pressure the maximum energy losses for protons, deuterons and 
tritons were calculated to be 0-48, 0:8 and 1-2 mev respectively; «-particles 
could lose their whole energy. The distributions of fig 2 show that the 
main peak at 76 cm pressure is due to protons. 


Counts 
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i ion 7Be(n p)’Li. The arrow 
height spectrum of protons from the reaction 

ie Pde the pulse Theight expected from protons leaving the residual 
7Li nucleus in the excited state at 478 kev. 


In fig. 3 the pulse height distribution is shown after the subtraction of 
background effects (e.g. particles from the 9B (n«) reaction fast Renpron 
recoils, radioactive contaminations etc.). Pulse height distributions 
were recorded for filling pressures of 76 cm and 11 em of mercury pressure, 
with and without cadmium in the incident beam, and ns ith the Be 
deposit facing the collecting volume and reversed. se — ae 
cedure yields a curve (fig. 3) which is the spectrum of pulses which can be 
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attributed to protons emitted from the 7Be as a result of slow neutron 
bombardment. 

A sample of 233U of the same physical dimensions could be substituted for 
the 7Be deposit, and the fission counting rate determined as a measure of 
the neutron flux. Six measurements of the proton pulse height distribution 
and the corresponding fission rate were made over a period of time. A 
proton counting rate was obtained by summing over the spectrum between 
1-1 and 1-8 mev and normalizing to constant fission rate, i.e. to constant 
incident neutron flux. This normalized rate is plotted on a logarithmic 
scale against time in fig. 4. The half-life of the decay of the proton rate 
is 50--2 days. The instability of the target material is thus established 


Normalized Fig. 4 
proton 
counting 
rate 


0 20 40 60 80 Days 


The decay of the 7Be(np) proton counting rate with time. 


and the value of the half-life is in fair agreement with that of 7Be. Further- 
more, if 10% of the initial proton rate were due to a background constant 
in time, the half-life would be reduced to 40 days and the fit to an ex- 
ponential decay worsened. It seems certain that the proton rate does 
not include such a background, i.e. the procedure for evaluating the 
background is unlikely to give too high a value for the (n p) cross section. 

Sufficient energy is available in the 7Be (n p) reaction for protons of 
reduced energy to be emitted and the residual 7Li nucleus left in the 
excited state at 478 kev. The arrow in fig. 3 indicates the position 
expected for the peak due to such ‘ excited state ’ protons. The proportion 
is less than 10%, of the ground state group. 
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In fig. 4, and the spectra from the other five runs, there is an asymmetry 
and a broadening of the peak, relative to the peak from the B (n «) 
reaction. The greater width of the proton distribution was shown to be 

due to the dependence of pulse height on collection time ; the contribution 
to the total proton counting rate of that portion of the spectrum below 
1-1 Mev was estimated to be less than 7°% from a study of the effect of 
collection time on pulse height distributions from «-particle sources. This 
method gives a better estimate than study of the pulse height distribution 
obtained by the subtraction procedure, where the statistical errors 
become large at the lower energies. 

The relation between the neutron flux and the fission counting rate of 
the ?33U sample was determined in a subsidiary experiment in terms of the 
capture cross section of gold (97--1 barns for neutrons of 0-025 ev energy). 
The accuracy of this measurement was +3%. 

The cross section for the (n p) reaction in *Be leading to the ground state 
of *Li, calculated from the proton counting rate per unit fission rate and 
the mass of 7Be is 5-30-8104 barns. Of the 15% error, 10% occur 
in the determination of the mass of 7Be because of the uncertainty in the 
branching ratio ; the other main source of error is in the interpretation 
of the proton spectrum. Ifthe proportion of small pulses has its maximum 
value, the proton rate may be 7%, too low ; on the other hand it may be 
10% too high if the proportion of protons leaving “Li in its excited state 
has its maximum possible value. 


2.3. The *Be (n «) Reaction 


As was discussed in the introduction, the parity of 7Be may be deduced 
from the relative magnitudes of the (n p) and (n «) reaction cross sections 
for neutrons incident with thermal energies. 

Some 40 hours were spent in examining the region of the energy 
spectrum between 9 and 10 mev with a gas filling pressure of three atmo- 
spheres. It is believed that traces of fissionable material were responsible 
for a background counting rate of 2 counts/hour/Mev energy interval. 
An upper limit of 1 barn can be assigned to the (n «) cross section on the 
assumption that a peak counting rate three times the standard deviation 
of the background rate would have been observed. 

The (np) reaction is therefore much more probable than the (n«) 
reaction, corresponding to odd parity for the ground state of “Be. The 
result is in accord with shell model predictions. 


§ 3. DESTRUCTIVE IRRADIATION OF 7Be 
3.1. Method 
In the presence of neutrons “Be nuclei are transformed by neutron 
capture as well as by K electron capture. The (n.p) cross section 18 so large 
that the increase of decay rate resulting from neutron bombardment gives 


a method of determining the neutron capture cross section. If the ratio 
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of the activities of two 7Be samples is measured to be y (> 1 say) and the 
stronger is irradiated for a time ¢ in a flux v, the ratio of activities after 
irradiation will be reduced to y!=y exp (—o fv dt), where o is the cross 
section for all processes in which 7Be is destroyed. For o=5 x 10-29 cm?, 
v=10!2 neutrons/em2/sec and t=2-6 10° sec (1 month) the change in 
relative activity y/y1—1—14%. The cross section can thus be measured 
in terms of a time integrated neutron flux and relative source strengths 
and does not depend on a knowledge of the 7“Be decay scheme. 


3.2. Results 


Two samples of 7Be were sealed off in quartz tubes and their y-activities 
compared with a scintillation detector. The value found for the ratio of 
activities (y) was 1-112-+-0-009. 

Fig. 5 
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Difference spectra from *7Be samples. x destruction of irradiated 7Be assumed 
to be zero. o destruction of irradiated *Be assumed to be 10%. Full 
line radiation from Co. 


The stronger of the two sources, together with a similar quartz capsule 
containing 10-37 mg of cobalt wire, was irradiated for a period of 28 days 
in a high flux position in B.E.P.O 

The irradiated *Be sample was found to have picked up a spurious 
activity of long life and y-ray energy ~900 kev. The spectra of the 
irradiated and unirradiated 7Be samples were studied with a sodium 


siete 
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iodide scintillation detector and thirty channel pulse amplitude analyser 
covering the energy region 300-660 kev. The value of the ratio of the 
"Be activities, y1, was determined as follows. The ordinates of the 

“spectrum from the unirradiated sample were multiplied by factors x 
before subtraction from the irradiated sample’s spectrum, to give ‘ differ- 
ence spectra’. The value of x for which the corresponding difference 
spectrum could most closely be fitted by the spectrum from a source of 
cobalt 60 (y-rays of 1-1 and 1-3 Mev) was taken as equal to y1. Identical 
values were obtained by fitting the difference spectrum to a **Mn spectrum 
(y-ray energy 850 kev) i.e. the value found for the ratio of 7Be activities 
was insensitive to the exact energy assigned to the impurity radiation. 
Difference spectra are shown in fig. 5. The upper set of points correspond 
to a 10% destruction of “Be by neutron irradiation, the lower set to no 
destruction. The full line is due to cobalt radiation. From four sets of 
data taken over a period of weeks a consistent value of y1=1-043-0-005 
was found. The same value was found when data were taken with an 
organic crystal detector. 

From the activity induced in the cobalt specimen the integrated flux 
( fv dt) was found to be 1-36--0-10 x 1018 neutron/cm?. 

Thus y=1-112-+0-009, y1=1-043-+40-005, fv dt=1-36-L0-10 x 1038. 

Hence o=4.8-+-0-9 x 10-?° cm?. 

If the spurious activity induced by neutron bombardment of the 7Be 
sample included a body of long half-life emitting y-rays not resolved 
from those of 7Be the destruction of 7Be will have been underestimated. 
If this possibility is admitted the quoted cross section must be regarded 
as a lower limit. 

With stronger 7Be sources and longer irradiation times the method is 
capable of considerably greater precision. The value found for the 
destruction cross section is in agreement with that measured for the 


(n p) reaction. 


§ 4. APPLICATION OF THE RECIPROCITY THEOREM 
4.1. Introduction 


To deduce the spin of 7Be from the reciprocity theorem (§ 1) it is 
necessary to determine the 7Be (n p) and 7Li (p n) cross section for energies 
of incident particles giving the same excitation of the compound nucleus 
in the two cases. This is not feasible experimentally ; a certain form of 
the energy dependence of the cross sections must be assumed. 

From the reciprocity theorem it is evident that near the threshold of 
the 7Li (pn) reaction, where the energy of the neutron 1s changing much 
more rapidly than that of the proton, a 1/v law dependence of the 7Be 
(n p) cross section is equivalent to a 7Li (p n) Cross section proportional 
to ./H,, where EH, is the energy of the neutron In the centre of mass 
system. Taschek and Hemmendinger (1948) show a curve of o/ JE, 
versus 4/H,, which falls steeply with increasing £,,. However, in the 
_ region H,,=17-30 kev, it appears that o/,/H,, has reached a constant value 
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of 0-045 barns/\/kev. These results are shown in the circles of fig. 6. 
On the assumption that o//H,, remains at this value as FH’, is decreased 
to thermal energies, the (n p) cross section for neutrons incident on *Be 
with an energy in laboratory co-ordinates of 0-025 ev is 1-6 x 10* barns, 
for a 7Be spin of 3/2, or twice this value for a spin of 4. The mean measured 
value of 5-0--0-8 104 barns corresponds to a spin of 0-14--0-15. It 
appears that o/\/H, cannot be constant. 


Go . 
barns/1/kev Fig. 6 
V (En) 


0 5 10 5 20 25 


E,, kev 
The *Li(p n) cross section near threshold. 
o points due to Taschek and Hemmendinger 
x points from *Be activation 
I points from proton recoil counting 


Three experiments have been made, with the proton beam of the Har- 
well electrostatic generator, in an attempt to clarify the situation: the 
absolute value of the 7Li (pn) cross section has been determined at two 
energies near threshold, the energy dependence of the (p n) cross section 
has been studied to within 5 kev of threshold and some observations made 
on the inelastic proton scattering. 


4.2. Measurements on the System *Li-+-p 


For the absolute cross section determinations the 7Be activity produced 
by bombardment with a measured charge of protons was compared with 
the activity of a source of Au standardized by coincidence counting. 
The superficial masses of the thin targets used were found from the neutron 
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yields far from threshold relative to the yields from targets thick enough 
to be weighed with accuracy. The values obtained for o/\/H, were 
5-24-0-6 x 10-6 em?/1/ kev and 4:8+0-5 x 10-26 em?/1/ kev for values of 
‘Lf, of 13 kev and 18 kev respectively. The major source of error arises 
from the uncertainty in the value of the branching ratio in 7Be decay ; 
the difference between the two values is probably significant. 

To follow the variation of cross section with energy rather nearer to 
the threshold the neutron yield from a thin target was measured with a 
hydrogen filled proportional counter set at 0° to the proton beam (Tunni- 
cliffe and Ward 1952). From the spectrum of proton recoils, the neutron 
energy in laboratory co-ordinates and the number of neutrons incident 
on the counter were determined. The neutron energy and the solid 
angle subtended by the counter in the centre of mass system were 
calculated therefrom. The total neutron yield was then derived on the 
assumption of isotropic emission of neutrons in the centre of mass system. 
The results are shown in fig. 6 and have been normalized to Taschek and 
Hemmendinger’s values at the high energy end. The errors increase as 
threshold is approached largely because the ratio of solid angles in labor- ’ 
atory and centre of mass co-ordinates is a rapidly varying function of 
energy near threshold; here also the shape of the proton recoil spectrum 
becomes extremely sensitive to target thickness and proton energy resolu- 
tion and the assignment of neutron energy and yield is untrustworthy. 

From fig. 6 the value of o/\/#,, is seen to be increasing as threshold is 
approached. The value of o/\/EZ, at threshold, determined from the 
7Be (n p) cross section on the assumption of a spin of 3/2, is 0-14 
barns// kev. 

The yield of 478 kev y-rays from inelastic proton scattering was studied 
near threshold. Below threshold there was no indication of a resonance ; 
above threshold the yield decreases due to competition from the (pn) 
reaction. It is noteworthy that the decrease was more abrupt than would 
be expected from a (pn) cross section proportional to \/H,,; a decrease 
of o/\/E,, with increasing #,, near threshold is also required to account 
for the behaviour of the inelastic scattering cross section. 


4.3. Discussion 


From the rather imprecise data on the ‘Li (pn) reaction it is not 
possible to obtain by extrapolation a trustworthy value of o/ \/H, as E,, 
tends to zero, i.e. the spin of 7Be cannot be determined from the reciprocity 
theorem. Some conclusions are possible if the procedure is reversed ; 
it is most likely that the mirror nuclei 7Be and “Li have the same spin 
(3/2) in which case o/\/H,, decreases by a factor 3 as LH, increases 
from 0 to 16 kev. If the variation of cross section is determined by a 
single level in *Be, Hy kev below the threshold and of width J’ kev it 
follows from the Breit—Wigner formula that 


3(H,2+I2/4)=[(Ey+ 16)?+ £7/4] 
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for s wave neutron emission. Hence J’ cannot exceed 30 kev, and the 
level must lie between 20 kev below the threshold and 2 kev above. The 
trend of o/4/Z,, with £,, suggests that the level is below threshold and of 
maximum width. Other levels may be contributing to the cross section, 
and neutron waves of higher angular momentum may be emitted at the 
higher energies. In the latter case neutron emission may not be isotropic 
and the derivation of the results of fig. 6 is invalidated. The effect of 
contribution from other levels, or from neutron waves other than s 
waves is that the value derived for the level width may be too large. 

The existence of a level in 8Be very close to the neutron threshold has 
been suspected for some time. Breit and Bloch (1948) suggested it in an 
analysis of the 7Li (p n) results of Taschek and Hemmendinger. Bashkin 
and Richards (1951) found a variation in the 7Li elastic proton scattering 
cross section which might be associated with such a level. More recently 
Newson and Gibbons (1954) have found more definite evidence for the 
level from a study of the “Li (pn) reaction. Their results, too, suggest 
that the level is below threshold. 

From the absolute magnitude of the (pn) cross section and the total 
width it is possible to calculate the neutron width. The reduced neutron 
width, calculated on the assumption that only s wave neutrons are 
involved, is 50 kev for a (1—) level and 90 kev for a (2—) level ; if the total 
width of the level is set equal to the proton width, the corresponding 
reduced proton width is found to be 20 kev. For a (pn) reaction between 
mirror nuclei the neutron and proton reduced widths are expected to be 
equal. Since systematic errors can arise in the extraction of reduced 
widths the agreement may be considered satisfactory. 

The reduced widths found are of order one hundred times less than 
single particle widths ; the level must involve a disturbance of the core, 
rather than result from the motion of an odd nucleon about the core. For 
a discussion of anomalously small reduced widths in the 1 p shell the 
reader is referred to a forthcoming publication by Lane and Thomas. 
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SUMMARY 


An event is described in which a K-meson decays at rest to a particle - 
at plateau ionization which has the properties of an electron of 
pBc~80 Mey. This brings to three the number of K-mesons in which 
the charged decay product is apparently an electron. 


$1. InrTRoDUCTION 


DvRING routine scanning of plates from Flight 30 of the Sardinia Flights 
of 1953, an event has been found which resembles those already described 
by Friedlander et al. (1954), and by Dahanayake et al. (1954). A heavy 
charged meson decays at rest. The charged secondary particle of 
ionization indistinguishable from the plateau value, had a visible range 
of 9mm, passing through three plates before leaving the stack. The 
mean value of the scattering parameter % was 0-27°/100u. This value 
will later be shown to be significantly different from that of a «-meson 
of the observed ionization. Accordingly, the event was identified as 
being a further example of the electron-decay of K-mesons, the first 
evidence for the existence of which has recently been put forward by 
Friedlander et al. (1954). The following are the criteria (Friedlander et 
al. 1954) which must be satisfied by a Ke event: 

(1) The estimate of mass of the primary particle must differ significantly 
from that of the 7-meson and also from that of the proton ; 

(2) the ionization in the track of the secondary particle must be 
indistinguishable from the plateau value (b*=1), while the estimate of 
pc must differ significantly from that of a u.-meson of plateau ionization ; 
and/or 

(3) if the track of the secondary particle is of sufficient length. there 
should occur a significant loss of energy by bremsstrahlung. In the 
present case no evidence was found for significant decrease of the energy 
of the secondary particle by bremsstrahlung, nor indeed would such be 
probable, since the visible length of trajectory is ~1/3 radiation length. 
For this reason, the identification of the secondary particle is based on 


the statistical argument (2). 


§ 2. PRIMARY PARTICLE 


The primary particle emerges from a star of type 10--4n. Four 
black tracks and one heavy fragment (Z~3) came to rest in the stack 
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without visible secondary particles. Four grey tracks other than the 
K-particle, left the stack. No attempt was made to trace the lightly 
ionizing particles. 

The primary particle passed through 3 plates before decay, having a 
mean range per plate 66404. From the estimated mean gap-length in 
the neighbourhood of the point of decay, it was concluded that the event 
occurred when the primary particle was at rest, or moving with negligible 
velocity. The total path traversed was 9-2 mm, corresponding, for 
mass 966 m,, to an energy of emission of 37-8 Mev, and a time of flight 
0-11 10-% sec. Mass estimates were made using the techniques (1), 
scattering-range and (2), ionization-range. 


2.1. Scattering-range 

The mass estimates were made by the constant sagitta technique 
using the scattering scheme D,=0-5y of Fay et al. (1954). They give 
the following relationship which holds to a good approximation for 
singly-charged particles 

M=273 x (0-5D,)* 

where WM is the value of the particle mass in terms of m,, D; is expressed 
in microns. The value of the exponent N given by Fay et al. is 2-32. 
We have, however, taken the value N=2-27 as being more appropriate, 
as it takes into account variations in the value of the scattering constant. 
whereas the value N=2-32 is based only on the assumed range—-energy 
relationship. Extensive calibration work to check these and similar 
tables using identified protons and z-mesons is at present nearing 
completion (Bristol—Dublin Group). 

Using the value V=2-27 the mass of the primary particle was found 
to be 1010772 770,: 

2.2. Mean Gap-Length—Range (@, R) 

The mass of the primary particle was also measured using the G/R 
technique (O’Ceallaigh 1953, 1954, Menon and O’Ceallaigh 1954). 

Two protons identified by 5-ray counts and one 7-meson were used as 
calibration particles. The observations are shown in fig. 1 on a log—log 
plot. The results have been fitted by straight lines by the method of 
least squares. The estimates of mass* are : 


From the 7-meson 789-+-100 m,. 
From the protons 867-++- 80 m,. 


. *Note added in proof.—A normalized @ calibration for the 8.30 stack using 
methods analogous to those developed by Fowler and Perkins (private com- 
munication) has since been completed. Values have been found for the 
normalized parameter G*=G)/@,. Gy is the plateau mean gap-length deter- 
mined by measurements on electrons, (, the mean gap-length of protons and 


a-mesons of selected range. Use of this calibration yields a mass-estimate of 
the parent K-meson of 970-+-70 m,. 
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These are in satisfactory agreement with those found for the range- 
scattering measurements. 
Thus, the decaying particle satisfies criterion (1) and is assumed to be 
‘a K-meson. 
Fig. 1 


MASS CALIBRATION 


MEAN GAP LENGTH 


RANGE (mms) 


§ 3. SECONDARY PARTICLE 


The trajectory of the secondary particle was flat and had a total 
length in 3 plates of 9mm. It left the stack through the processed edge 
of the fourth plate. 

Since the plateau value of blob-density was very low (14 blobs/100,.), 
and the significance of the results depends on the accuracy of the follow 
through procedure, a very precise technique was necessary. This was 
similar to that at present employed in Bristol (1954). A C.T.S. M.4000 
microscope was employed, having a rotating turntable carrying an 
accurately constructed jig secured by brass clamps. By use of this, the 
position of the track was fixed to within a few microns. It was orientated 
along the x-axis of the stage and the point of entry in the succeeding plate 
was calculated from the dip of the track. The spurious shift along the 
y-axis was found by examining under low power the displacement of 
steep neighbouring tracks. The stage was set in the calculated position, 
and the continuation of the track located. The identification. was 
verified by comparison of the track characteristics, 6, « and dip angle. 
This operation was checked by several observers. Owing to the care 
taken, it is felt that the possibility of a mistaken identification is excluded. ° 


3.1. Influence of Distortion 
The scattering measurements were continued up to 1mm from the 
edge of the third plate. It was necessary, therefore, despite the small 
dip of the track, to estimate the influence of distortion on the measured 


2H2 
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values of scattering parameter. In this plate, 7(a), the surface distortion 
vector component at right angles to the track was ~30, corresponding 
for 150 cell, to a spurious sagitta of 0-005 on a track of length 3500p in 
a 600 plate. This would lead to an overestimate of ~1% of the « in this 
plate. In the other two plates, 6(a) and 5(a), both the distortion vector 
and dip are rather less than in plate 7(a). Accordingly, we may neglect 
the influence of distortion on the estimate of p8c of the secondary particle. 

Stage noise has been measured for the instrument used and shown to 
be negligible for cell sizes <500p. 

The results of the scattering and ionization measurements on the 
secondary particle are given in the following tables. The measurements 
have been made with coefficient of overlap A=2 and 3 but NV throughout 
refers to the number of non-overlapping cells in each sample. 

Table 1 contains details of measurements and characteristics of the 
individual segments of track. 


7(a) 3536 12-6+0-6-| 13-29+0-27 0-95 


+ 0-05 


150/50] 0-304 | 20 
+0-055 


Table 1 
Projected Plateau 4 Cell a 100 . 
Plate | fength (u) | POPs OOH | Piohsi100 pn) 2 |sizeata)| (OOO mat 
5(a) 1820 14-4+-0-7 | 13-82-+0-21 1-04 | 75/25 OS Tise20 
+0-05 +0-056 
6(a) 12:-8+0-6 | 13-19 +0-26 0-97 | 100/25 0-300 | 33 
+0-05 +0-042 


Results of more extended measurements on the complete track using 
the ordinary method of noise-elimination (Fowler 1950) are given in 
table 2 

The weighted mean values of pfce were as follows : 


81+ 9 Mev (a(2)), 
ppc 98-+-11 Mev (a), 


Appropriate values of scattering constant (Gottstein e¢ al. 1951, fig. 1) 
have been used to calculate the values of pc. 

It will be observed that method of third difference yields consistently 
lower values of the estimate of a5). It would appear doubtful, therefore, 
if an unbiassed estimate of second difference results from the division of 
the mean value of the third difference by the accepted factor 21/2, 
Evidence obtained from extensive independent measurements, indicates 
that this factor should be decreased by ~7°%, (Appendix). It appears 
reasonable, therefore, to decrease the estimate of pBe from «®) to 
90+ 11 Mev. 


2 aes eon amends 
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renee eS 


(A) Rigorous 4 D cut-off 


ee a eS a ee 


Cellsizes(z)| 150-100 | N 150—50 N 100—50 N 
a2) 0-311 +0-035 52 | 0:324+0-035 | 52 | 0:268-+0-024 | 80 
pc (Mev) 76+8 73 +8 86+8 
a3) 0-243 +0-027 48 | 0-246+-0-027 | 48 | 0-232+0-022 | 76 
Bc (Mev) 98 + Ja 96 SI 99 +9 


(B) 4 D cut-off with replacement 


at?) 0-325 +0-036 | 52 | 0-314+40-035 | 52 | 0-298-+0-027 | 81 
pBe(mev) | 75-549 780-49 80:0-£7 

al) 0-249 +0-029 | 49 | 0-249+0-029 | 49 | 0-260+40-023 | 81 
pc (Mev) 98411 he sO8 2277 9248 


(C) No Cut-off 


5) 0-313 40-035 | 52 | 0-32040-036 | 52 | 0-294+0-027 | 81 
pBc (Mev) 83 +49 S1ae9 86-5 +8 

a3) 0-248 40-028 | 49 | 0-245-.0-247 | 49 | 0-253+40-023 | 78 
pBc (Mev) 105-412 106 +12 100-+9 


3.2. Calibration b*—pBc 

In earlier work it has been customary to rely on the results of Daniel 
and Perkins (1954), Daniel et al. (1952) which give the relationship 
between g* and pfc for identified 7-mesons and protons emergent from 
high-energy interactions. This work was carried out by simultaneous 
blob-counts and scattering measurements. The ‘ best z-line’ of Daniel 
and Perkins was taken as standard, and it has been used to identify fast 
particles secondary to K-decay. 

Recent work at Bristol has indicated, however, that the values of g* 
for a given pfc may differ significantly from stack to stack. Consequently, 
it has been felt necessary to establish experimentally this relationship for 
the present stack. This calibration has been carried out by making 
blob-counts on the tracks of long flat 7-mesons which underwent 7-p-e 
decay or gave rise to characteristic o-stars. Five such m-mesons were 
measured in the interval of residual range 0-8 cm corresponding to pBc 
0-120 Mev. The relationship between Rk and pfc was taken to be that 
of Fay et al. (1954), points beyond the range 5-2cm being found by 
graphical extrapolation. For each plate the blob-count corresponding 
to plateau was found by observation on long flat electrons from. 7—p~-€ 
decay. This information was supplemented by counts” on similar 
electrons from p-e decay, which were accepted when they did not differ 
significantly from each other, and from those from 7--e decays in the 


NORMALISED BLOB-COUNT 
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same plate. Observations were confined to those portions of the tracks | 
which lay within 30. from the surface and 10 from the glass. Approxi- 
mately 1500 blobs corresponding to 1 cm of track were counted in each 
plate. No attempt could be made to assess variations of by with position 
in the plate, except that measurements were accepted only if made on 
segments of track more than 2 mm from the processed edges of the plate. 

The results are given in fig. 2. They are well fitted for 
120>pfc>20 Mev by a straight line, log b+=1-260—0-5626 log pe. 
The extrapolation of this line to cut the line b*=1. yields an under- 
estimate of the value of pc at which a z-meson first reaches the plateau 
blob-count. The line agrees with the ‘ best line * of Daniel and Perkins, 
at the point pSc=100 Mev, but appears to have a steeper slope. The 
p-line was drawn in using the mass ratio m,/m,—=1-:321. The points of 
intersection of these lines with the line b*=1 were found to be ] 

: 
| 
| 


pBc=173 Mev, 
pBc=131 Mev. 


Fig. 2 
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§ 4. CONCLUSIONS 


As shown the value of pe for the secondary of the present event was 
81+-9 Mev and is 6 standard deviations from the expected point. Thus, 
the probability is negligible that the secondary particle is a p-meson. 
Even accepting the results of the estimate of pfc by the method of third 
differences, the experimental pc is significantly different from that 
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corresponding to ay-meson. We may therefore conclude that the present 
' example is to be identified as a case of K-e decay bringing the number of 
those identified at present to 3. 

The observed value of pc agrees well with that found by Friedlander 
et al. but differs significantly from that found by Dahanayake et al. 
(49+-4 Mev). There is therefore good reason to suppose that in addition 
to the electron, at least two neutral particles are emitted as products of 
decay. This conclusion would be in agreement with the predictions of 
Michel (1952) based on the hypothesis of a universal Fermi interaction, 
analogous to the fundamental process of u-e decay. 

An event which from preliminary measurements appeared to be a 
second example of a K-e process has since been very carefully remeasured. 
The primary K-particle of estimated mass by scattering-range 970-L!%? 
decays at rest to a secondary particle of blob-density indistinguishable 
from plateau. This was traced through 11 plates having average pro- 
jected length 900, per plate. When however distortion corrections are 
applied, the estimate of pfc of this particle becomes 120+8 Mev. Thus. 
it appears as likely to be a »-meson as an electron, although it is most 
improbable that it is a 7-meson, and hence the reference (Dahanayake 
et al. 1954) requires modification. 
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AP PEND TeX 


It has usually been assumed that an unbiassed estimate of D, may be 
obtained from the experimental value of D, by dividing by the factor 2", 
This belief is founded on the assumption that the difference between two 
independent normally distributed variates is itself normally distributed. 
with variance equal to the sum of the parent variances. However, the 
individual values of D, taken even with coefficient of overlap A=1 are 
not independent, and also their distribution is only approximately normal. 

Our measurements by the Cosyns-Vanderhaege method showed that 
the distortion of the secondary of the K-particle described in this paper 
was negligible. We have been led to suspect therefore, ED the source 
of the discrepancy between the estimates of «199 obtained from measure- 
ments of second and from third differences, might be a bias introduced by 
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non-satisfaction of the conditions of normality and independence. Hence, 
the accepted value 21/2 may be inappropriate to the experimental 
circumstances. 

For these reasons we have investigated the matter experimentally 
using three groups of tracks with the results shown in table 3. 


Table 3 
Nature of tracks No of cells D,/D» 
m-mesons (constant sagitta) 1806 1-26 40-034 
Protons (constant sagitta) 612 1-25 +0-064 
105 Mev positrons from Cornell synchrotron 1449 1-35 +0-040 


From table 3 it would appear that the appropriate value of the divisor 
is lower than 21/2, and we have adopted provisionally the value 1-32, 
corresponding to 7°, bias. Investigation of this effect is being continued, 
because of the apparent discrepancy between the results for constant cell 
and constant sagitta. 
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SUMMARY 


Data on Brillouin Zones for nine aluminium-rich compounds of iron 
have been obtained by x-ray analysis. For binary compounds and a 
ternary compound with nickel, the data can be correlated by assuming 
a reasonable electron absorption by iron, in agreement with previous 
results for alloys of cobalt and nickel. For ternary compounds containing 
copper and silicon, such correlation is not possible as the absorption 
required is obviously too high. Similar difficulties are encountered with 
some manganese compounds. It is suggested that the presence of a set 
of very strong x-ray reflexions at nearly the same spacings 1s not a 
sufficient criterion for classifying these complex phases as electron 
compounds. 


§ 1. INTRODUCTION 


THis paper is concerned with the application of Brilloum Zone analysis 
to aluminium-rich intermetallic compounds of iron. Many intermetallic 
phases have been classified as electron compounds because it is believed 
that their stability is critically dependent on the filling of the zone by the 
free electron states. Brillouin Zone planes are related to strong X-ray 
reflexions (Jones 1934), and this relationship has been used to calculate 
volumes of Brillouin Zones for aluminium-rich compounds of the metals 
from chromium to nickel. It has been assumed that these are electron 
compounds and the electron—atom ratios derived from the volumes of 
the Brillouin Zones have been combined with the compositions of the 
compounds to give values for the ‘valency’ of the transition metal 
atoms. This analysis has been applied by Raynor and Waldron (1949), 
Douglas (1950), Robinson (1952 a, b, 1953), Robinson and Black (1953) 
and Black (1954) and these results and other work on these compounds have 
been reviewed by Taylor (1954). In general, valencies in approximate 
agreement with the electron absorption hypothesis of Raynor (1944) have 
been obtained. In the present work an attempt has been made to review 
an extensive set of such data for iron compounds. 


= Waerlheee 
* Communicated by W. H. Taylor. 
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§ 2. ANALYSIS OF THE DATA 


It is necessary to consider precisely the significance of the data and this 
can best be done by reference to the work of Jones (1937) on the «/B phase 
boundaries, in which he analysed the role of the Brillouin Zone in determin- 
ing the relative stability of these two phases. Curves of N(#)d# against 
E (where E is the kinetic energy of the electrons) were computed and 
from them the electron—atom ratio at which the kinetic energy of free 
electrons is the same for both structures was found; on either side of 
this value, one or other of the two structures would be the more stable 
because the electrons could be contained in the metal with a smaller 
total energy. The zones for the two structures are not good approxima- 
tions to a sphere and their inscribed spheres contain only 64%, and 78%, 
of the full zone volumes. Jones’ results show that each N(/) curve falls 
sharply after the Fermi sphere first touches a plane of the zone, so that 
although both zones have a full volume corresponding to a ratio of two 
electrons per atom, the f-structure becomes more stable than the « 
when the latter has an electron—atom ratio of 1-41. In the present work, 
the largest inscribed spheres contain about 80°, of the zone volume. 
The N(£)dE curve will begin to fall after the Fermi surface first touches 
a plane of the zone, and some alternative structure should thus be more - 
stable at an electron—atom ratio just above that required to fill the 
inscribed sphere. This latter value will be quoted here, together 
with the electron—atom ratio for a full zone—which should be an upper 
limit, for beyond this the electron distribution must overlap the zone 
faces and give a large increase in the energy of the electrons. 

In most of the results given below, zone volumes have been calculated 
from the areas of the zone faces as determined graphically from drawings 
of the zone in projection. In other cases, an estimate of this volume 
has been made, and these figures are quoted in brackets. To each value 
of the electron—atom ratio derived from the zone, there corresponds a 
value for the number of electrons that each transition atom must absorb 
if the compound is to have this ratio and if normal valencies are assumed 
for aluminium, copper and silicon. 

In all cases, the selected reflexions form an outstandingly strong group 
of approximately constant p-value. The p-values either are tabulated 
here or can be found in the original sources. In tables 1, 2 and 3, an 
approximate estimate is given, for each compound, of the ratio of the 
mean intensity of zone reflexions to the intensity of the strongest reflexion 
which appears at lower angles, i.e. inside the zone. In all cases, there 
are only one or two inner reflexions which are not very weak, and these 
do not form an approximately spherical zone. The figures may be 
compared with results for y-brass, a complex structure for which it is 
generally accepted that the zone defines a stable electron concentration. 
In this case, the {222} reflexions lie within the zone formed by the {330} 
and {411} reflexions, and the ratios of the intensity of {222} to the 
intensities of {330} and {411} are approximately 1:9 and 1 : 3-5 
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respectively (Bradley and Thewlis 1926).* The zone reflexions themselves 
have unitary structure factors (Lipson and Cochran 1953, p. 251) of 0-7 
and 0-3 for y-brass: for FeAl,, FeCu,Al,, Fe,Al, and 8(AlFeSi) for 
which these unitary factors can be calculated, the values lie in the range 
0-5 to 0-7. It seems that any criterion for restricting the application of 
zone analysis which is based on the relative or absolute intensity of zone 
reflexions cannot distinguish between these compounds and y-brass. 


§3 
3.1. The (Al-Fe-Si) System 
The existence of three ‘«-phase’ compounds has been reported 
(Armand 1952) and results for two of these have been discussed by 
Robinson and Black (1953). Chemical analyses by different workers of 


extracted crystals of the B-phase have given the following results :— 
Fe(wt.%)  Si(wt.%) 
Dix and Heath (1928) 27 15 
Phragmen (1950) 27:3 13-5 
Pratt and Raynor (1951) 26-7—27-3 13-8-14-9 

The best approximate composition is Fe,Si,Al,, which corresponds to 
27-2 wt.% iron and 13-7 wt.% silicon. 

Phragmen reported that the phase is monoclinic with a=b=6-124, 
c=41:54 and B=91°. For the present work, crystals obtained by Pratt 
and Raynor were examined by Laue, oscillation and Weissenberg 
photographs, and were found to possess Laue symmetry 4/m, with the 
following cell dimensions :— 

a=6:18+0-06A4 u. c=42-5-+0-54 u. 
The Brillouin Zone was obtained from the following strong low-angle 


reflexions :— 


Indices | Multiplicity | p=1/2doay (A) 


h 

Dea ee 0) 4 0-229 
1 tk ae 8 0-241 
Omer 1.6 8 0-248 
ioe - 8 0-257 


The density, measured by a flotation method, is 3°50-+-0-10 g/cm?, 
which indicates 8-3 units of Fe,Si,Al, per cell. The Brillouin Zone 
data are summarized in table 1. The mean atomic volume is obtained 
from the measured density and chemical composition. 

A third phase, the 5-phase, has been observed in this system by Gwver 
and Phillips (1937), Janiche (1936) and Phragmen (1950). Phragmen 


= mm syylic ic 
* Relative intensities are considered here, whereas Bradley and Thewlis list 


structure amplitudes. 
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calls it the t-phase ; he has measured the unit cell, which is tetragonal 
with a=6-124 and c=9-484. According to Janiche, the composition is 
Al,Si,Fe. An approximate structure has been obtained by Spiegelberg 
(1954), which indicates the composition Al,Si,Fe, and from his results it 
has been calculated that the following low-angle planes give very strong 
reflexions :— 


: Indices Multiplicity | p=1/2d¢xn (A-4) 
eet. 4 0-231 
1 1 4 8 0-241 
ee es sf 0-258 
ee eee, 4 0-267 


The mean atomic volume is obtained from the structural data and unit 
cell dimensions. 

The data for this system are summarized in table 1—results for the 
binary compound FeAl, are also included (Black 1954). The «,-«.-x, 
series might be considered in terms of an increase in silicon content being 
balanced by an increase in transition metal content to maintain a constant 
electron—atom ratio. The 8 and 6 phases do not fit this trend and would 
have higher electron—atom ratios. This increase is not reflected in ap 
increase of Brillouin Zone volume. 


3.2. The (Al-Fe—Cu) System 

The compounds occurring in this system have been investigated by 
Phragmen (1950): there are two ternary compounds. Phragmen 
reported that the first of these is isomorphous with MnAl, ; its composi- 
tion, in atomic percentages, is 83-2°, aluminium, 4:2°, copper and 
12-5%, iron, compared with 14:3°, manganese and 85-7°, aluminium in 
MnAl,. A binary compound FeAl, has never been found and it is difficult 
to explain the role of copper in stabilizing a compound with a composition 
approximating to this. If the cell has the same size and symmetry, 
then there is not sufficient copper for complete replacement of aluminium 
or manganese on one particular site in the unit cell. Phillips (1954) has 
suggested that a compound FeAl, cannot be formed because the short 
Mn-—Al] distances of MnAl, are not possible for a transition metal which 
does not absorb as many electrons as manganese and that copper relieves 
this strain by replacing the aluminium. In FeAl,, however, (Black 1954) 
the Fe—Al distances lie in the same range and can be as short as the 
Mn—Al distances of MnAl,. The Brillouin Zone and corresponding 
electron—atom ratios for the compound can be calculated from the 
results for MnAl, (Nicol 1953) and the cell dimensions given by Phragmen ; 
his measured density, composition and cell size indicate that the mean 
atomic volume is almost the same as in MnAl,. The data are given in 
table 2. Several compounds of this kind are related in such a way that 
isomorphous compounds have the same electron—atom ratio (6(AIMnSi) 
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and Co,Al;—Robinson (1952 a), FeNiAl, and Co,Al,—Bradley and 
Taylor (1940)), but this type of explanation cannot apply to this case. 
Only if iron has about the same effective valency as manganese can 
_ these two compounds have the same electron—atom ratio. 

The second ternary is of composition FeCu,Al, and has been examined 
by Wiehr (1940) and Westgren (see Phragmen 1950). The unit cell is 
tetragonal with a=6-32A, c=14-78A. The structure was analysed by 
Westgren and the parameters have been published by Phragmen ; from 
these the structure factors for all planes with sin 6/A<0-5 have been 
calculated. The largest of these are listed below :— 


, indices | | Multiplicity | p=1/2daan (A) 
Demet d 16 0-223 
oe. 2s 16 0-245 
Beet FO 8 0-250 
Sey 16 0-253 
0 0 8 2 0-271 


| 
The mean atomic volume is obtained from the structural composition and 
the cell parameters. 

The data for phases occurring in this system are reviewed in table 2. 
The trend of composition shows that the iron content is reduced by 
replacing both iron and aluminium by copper and this is consistent with 
a tendency to maintain a constant electron—atom ratio. CuAl, has been 
included for purposes of comparison. This structure gives strong low 
angle reflexions with values of p ranging from 0-211 to 0-263. It could 
not be classified as an electron compound—the normal Cu and Al valencies 
indicate that the electron: distribution must overlap the Brillouin Zone 
planes. 


3.3. Al-Transition Metal Compounds 
Data for FeAl, have been presented elsewhere (Black i954) and are 
quoted in table 1. The structure of another binary compound, Fe,Al, 
has been determined by Schubert (1953), and Brillouin Zone data have 
been calculated from his results. The unit cell is orthorhombic, with 
a=7-68A, b=6-40A and c=4-204 and space group Dij-Cmem. The 


strong planes are as follows :— 
SSS 


‘ me Multiplicity | p=1/2dcan (A“) 
Soon ih 8 0-235 

ee he? 2 0-237 

oy geal 8 0-241 | 
ae = 0 4 0-243 

ih eal ae 8 0-258 

4 0 0) 2? 4-260 
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Schubert’s result indicates that the composition Fe,Al, is attained 
by having some of the atomic positions only partially filled. The mean 
atomic volume is calculated from the structural data and so makes 
allowance for the effect of this. Bradley and Taylor (1940) have reported 
the existence of a compound Fe,NiAl,), which is isomorphous with 
Co,Al, ; from calculations of the zone for this latter structure by Robinson 
(1952 a), and assuming the same cell parameters for the two compounds 
(Bradley and Taylor reported only small differences), the results of table 3 
have been obtained. In estimating the iron absorption, it has been 
assumed that nickel absorbs two electrons less than iron, and the com- 
position at the electron-rich limit of the Fe,NiAl,, phase field has been 
taken. 


§ 4. CONCLUSION 


If it is assumed that a structure has all interatomic distances approxi- 
mately equal, and high co-ordination so that each atom is completely 
surrounded by neighbours, then it can be shown that it must give a set 
of strong reflexion planes. The interplanar spacing of these can be 
related to the interatomic distance, by using the transformation which 
relates the intensity distribution in a powder pattern to the mean radial 
Patterson function (see formula (150) of Cochran and Lipson (1953)). 
It can then be calculated that for any structure of this type, the limiting 
Fermi sphere will contain. 1-5—2-5 states per atom. Only general atomic 
packing conditions are assumed in deriving this, and the result is inde- 
pendent of the scale of the structure as the product (mean atomic 

_volume < Brillouin Zone volume) is dimensionless. If empirical evidence 
from Brillouin Zone correlations is to be significant, it must show that the 
individual deviations within this mean range correspond to differences in 
the compositions of the compounds. In many cases, a detailed correlation 
of this type is not found. 

Most of the results in 3.1 and 3.2 show that the enpirical correlation 
of zone volumes on the basis of a reasonable assumed valency for the 
transition atoms is not possible ; the results of 3.3 might fit approximately 
with the Raynor hypothesis and in this respect they are in agreement 
with results for other transition metal-aluminium compounds. There 
are not enough data on the ternary compounds with copper or silicon of 
other transition metals to assess whether the discrepancies found here 
are peculiar to iron. Reference may be made to table | of Taylor (1954) 
for previous results : the data on an aluminium—manganese—zinc ternary 
compound and on the binary (MnAl,) compound would both indicate 
a very large absorption for manganese. Of the two ternary compounds 
of manganese with silicon, the zone for the B-phase indicates an absorption 
of about 3 electrons for manganese, but the «-phase results again require 
a very much larger absorption (Robinson 1951). Jones (1953) has 
suggested that absorption to the extent of four electrons may be possible 
if this is compensated by a decrease in the free electron density im the 


neighbourhood of the atom. 
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It might be suggested that some of the compounds examined above do 
not belong to the same class as the aluminium-rich binary compounds. 
An attempt must then be made to define a criterion for distinguishing 
such compounds. The results given here suggest that ternary compounds 
containing silicon and copper do not belong to the same class as the 
binary compounds. An examination of their structures shows, however, 
that a distinction of this sort is not indicated by the inter-atomic distances 
and co-ordinations. It can be shown that the structures of several of 
the compounds discussed above are related, and that, in common with 
other compounds of manganese, cobalt and nickel, they all possess 
certain characteristic features. This will be discussed elsewhere (Black 
1954). 

The theoretical justification for extending the use of the Brillouin 
Zone theory to these complex compounds has been examined critically 
by Vousden (1954). At the same time, it seems that it is more reasonable 
to account for their stability in terms of a strong interaction between 
aluminium and transition metal atoms (e.g. FeAl,, Black 1954), and 
Vousden has shown how this approach explains their characteristic 
properties (1954). 
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SUMMARY 


Some results previously reported by the author are examined on the 
assumption that the thermal acceptors are copper atoms. A thermo- 
dynamic treatment of the solid solubility of copper in germanium is 
given in terms of the Principle of Lowest Free Energy, enabling a physical 
picture of the solubility to be given. This explains the approximately 
exponential dependence on inverse absolute temperature (below 800°C), 
which formerly led to the belief that thermal acceptors were lattice 
defects. An equation is derived which relates the concentrations of 
copper on the liquidus and solidus to the temperature. The experimental 
results are shown to obey this equation, which is equivalent to the one 
previously derived by Thurmond and Struthers, using a different approach. 


§1. IyTRoDUCTION 


Iv has long been known (Torrey and Whitmer 1946) that if a germanium 
specimen is heat treated at a high temperature, and then quenched 
rapidly to room temperature, the concentration of acceptor centres in it 
is increased. The concentration of these added acceptors, which are 
called. ‘thermal acceptors’ is determined by the temperature of heat 
treatment. 

There have been several publications (Fuller and Struthers 1952, 
Slichter and Kolb 1952, Seiler, Geist, Keller and Blank 1953) which 
suggest that the formation of thermal acceptors is associated with the 
diffusion of impurities into the germanium during heat treatment. 
In particular, copper has been suggested as being responsible for the 
observed heat treatment effects. Fuller and Struthers (1952) published 
measurements of the saturation concentration of diffused copper in 
germanium heat treated at various temperatures, and of the concentration 
of thermal acceptors in the same specimens. They found that within 
the limits of experimental error, the number of thermal acceptors was 
equal to the number of copper atoms. They therefore interpreted their 
results to mean that thermal acceptors were to be identified with copper 
atoms, and that the equilibrium concentrations of thermal acceptors 
were the saturation solid solubilities of copper in germanium at the 
temperatures in question. They found that the concentrations of thermal 
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acceptors and of copper atoms were given by relations of the form (1) 
apart from those at the highest temperature of measurement (919°C) 
which were lower than this relation would predict. 


n=A exp (— zp) eee ne Aeererg ( L) 


where n=concentration of thermal acceptors or copper atoms, 
A=constant, 
k=Boltzmann’s constant, 
T=absolute temperature, 
W=constant (=28 000 cal/mol. or 1-24 ev/atom). 


> 


7 


§2. EXPERIMENTAL WorK 


Measurements of the concentrations of thermal acceptors have been 
carried out in this laboratory (Hodgkinson 1954), and have been extended 
up to 930°c. 

Specimens of germanium about 0-1 x 0-12 cm? were heat treated in 
evacuated silica capsules placed inside an electrical resistance furnace. 
At the end of a period of heat treatment, the silica capsule was pulled 
from the furnace into a beaker of water. The total time of quenching 
was of the order of a second or less. Concentrations of thermal acceptors 
were determined from electrical resistance measurements, carried out 
in a thermostat at 30°c. 


2.1. Results 


The results of the measurements are shown in fig. 1, where Inn is 
plotted against 1/7. In the temperature range below 800°c, the con- 
centrations of thermal acceptors are given by a relation of the form of 
eqn. (1). The values of W determined for the individual specimens 
vary from about 1-13 to 1:57ev. Utilizing all of the measurements 
reported in this paper, n=1-61x 10% exp (—1-31/kT). Above 800°C, 
the concentrations of thermal acceptors are lower than those given by 
eqn. (1), and above 880°c, the concentration of thermal acceptors decreases 
with rising temperature. The measurements at 919°c made by Fuller 
and Struthers are consistent with this decrease. 


§3. THEORY 


If thermal acceptors are copper atoms, it would be expected that their 
equilibrium concentration would be determined by equilibrium in the 
system copper-germanium. Figure 2 shows the copper-germanium 
phase diagram (Smithells 1949). The thermal acceptors would form 
a solid solution of copper in germanium at the right hand side of the 
diagram, but the saturation solubility is too small to be shown. It will 
be seen that above 650°c, the solid solution will be in equilibrium with the 


liquid phase. It would be expected that the solid solubility of copper 


2E2 
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in germanium would fall to zero as the melting point of pure germanium 
is approached. However, it is the retrograde solubility of copper in 
germanium at lower temperatures which is of main interest here. 


Fig. 1 


Thermal acceptors/cm? 


T 


Experimental values of the variation of thermal acceptors with 
temperature, 


3.1. Previous Theoretical Treatment 


Recently, Thurmond and Struthers (1953) have treated the problem 
theoretically, by equating the partial free energy of copper in solid 
solution in germanium to that of copper in the liquid Cu-Ge phase in 
equilibrium with the solid. An expression is derived for the variation 
with temperature of K, the partition coefficient of copper between solid 
and liquid germanium. 
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4HJ—4H; , 0 AS¥ 


In K= RT roe aes a In Vor ; ° 5 (2) 


where 


K=ratio of concentration of copper on the solidus to concentration of 
copper on liquidus, 

AH,'=heat of fusion of copper (assumed independent of temperature), 

4H,*=partial molar heat of solution of solid copper in solid germanium, 


o=entropy term to account for changes of vibrational entropy of 


copper atoms on being transferred from a copper to a germanium 
crystal, 


AS,'=entropy of fusion of copper (assumed independent of temperature), 
yol=activity coefficient of copper in the liquid Cu-Ge phase, 
R=gas constant/mole. 


Fig. 2 
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The solid is treated as a regular solution, and initially the liquid is 
treated as ideal. Later, the activity coefficient y,’ is introduced to take 
account of the departures from ideality of the liquid, and the term o 
is introduced to take account of the change in vibrational energy of the 
copper atoms on being transferred from one medium to another. 
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Thurmond and Struthers conclude that In K should be a smooth 
function of temperature. On plotting experimental values of In K 
against 1/7’, a straight line is obtained. Extrapolation to 7=melting 
point of pure germanium yields a value of K=0-98 x 10~°, which is in 
good agreement with the value 1510-5, found by other methods. 
On the assumption that In y,! has little effect on the slope of In K v. 1/T, 
the slope is equated to 4H,’—AH,’, from which 4H,° is calculated 
(=45 200 cal/g atom) since 4H, is known. 


3.1.1. Present Treatment 

The treatment presented here differs from that of Thurmond and 
Struthers, by considering the integral free energies of the solid and liquid 
phases. The principle of lowest free energy (Lipson and Wilson 1940, 
Cottrell 1948) is applied to obtain the equilibrium conditions. This 
enables a simple physical picture to be built up, which gives a qualitative 
explanation of the form of the solid solubility curve of copper in 
germanium. To apply the principle of lowest free energy, the variation 
of free energy with composition for the two phases is calculated. The 
theory is shown to be consistent with the experimental results. 


3.2. Principle of Lowest Free Energy 
Equilibrium conditions in an alloy system may be considered in terms 
of the principle of lowest free energy. To apply this principle the form 
and relative position of the free energy curves of the two phases in 
equilibrium must be known. These are now investigated for the solid 
solution of copper in germanium, and for the liquid copper-germanium 
phase. 
3.3. Free Energy Curve for the Solid Solution 
It will be assumed that the solid solution of copper in germanium is 
substitutional ; however, the treatment for an interstitial solid solution 
will be similar, and the expression for the free energy will be of the same 
form. 
The free energy, f, per atom of a dilute substitutional solid solution 
is given by an expression of the form (Mott and Jones 1936) 
; kT N! eG dL: 
lene N A (N—n)!n! pe Mi ree) 


where N—=number of lattice sites/cm’, 
n=number of solute atoms/cm’, 


Frome 20) 


“0 0 


7’=absolute temperature, 
C',=heat capacity per atom at constant volume, 
H(n)=internal energy per atom at absolute zero of temperature. 


The first term k/N In N !/(N—n)!n! is the configurational entropy of 
mixing of the solvent atoms and the solute atoms. 

For a very dilute solution, each solute atom will be surrounded com- 
pletely by solvent atoms, and the energy increase of the solid due to the 
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addition of one solute atom will be independent of n. E(n) may therefore 
be written 


nH 
(n)=By+ =, 


where Hy is the value of H(n) for the pure solute, and £ is the energy 
increase oi) the solution on introducing one solvent atom, substitutionally. 
In addition, in the case considered here, the solute atoms will be ionized, 
and an additional term will have to be added to the free energy (Goldberg 
1952) in eqn. (3). Reiss (1953) has dealt with the problem of impurity 
atoms in semiconductors, and has shown that they may be treated as an 
ionized electrolyte. 

At the temperatures considered here, the germanium will be intrinsic 
and the impurity atoms will all be ionized. The extra energy term will 
thus be independent of temperature. Suppose it is equal to «, and 
suppose H+e=W. In addition, an entropy term k7'In2 must be 
added to account for the two possible spins of the electron in the impurity 
levels. 

If it is assumed that germanium is an Hinstein type solid, then for the 
pure solid (Mott anu. Jones 1936) 

== ——— In, 3 « Gy. . (A) 


i NG eel 
where v is the vibration epee * the atoms. 
If n atoms of solute replace n germanium atoms, they will have 
a different vibration frequency v’ from the atoms they replace. Suppose 
that the vibration frequency of the x nearest neighbours to a solute atom 
is changed to v” along a line joining them to the solute atoms. Then 
eqn. (“) becomes (Mott and Gurney 1948) 


£O oo 3kT , hv 


tC a = 
nl {3N —n(a+3) }kT es a 
hy u" 
+3nkT In 7 tna In or |: ws 50(5) 
Thus, eqn. (3) becomes 
! Wn nkT 
i meninia ee NO 


i hv hy" 
ef WV | oN mets yyer In ai +3nkT In 7 bnaklD In or | eG) 


Applying Stirling’s approximation to the first term on the right-hand 
side of eqn. (6), and aaa with respect to n, 


hy hy’ 
N 2 =e! In + > aa W—(a+3)kT In PP —kT \In 2+-3kT In — ip 


“of Vine 
—— = W—kT InO—kTin2. . {7) 
ok In TT =kT ln ———> a + 


where C= jayne 
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. ofa 
For a minimum, an ==0 
W W 
L.€: n=20(N—n) exp (- rr) ==2C'N exp (- ir) : 
Suppose n=n, at the minimum, then 
W 
Np=—A exp (- ir) Rie a kt 
where A=2C(N—n)=2CN= const. 


To calculate the solid solubility curve of copper it is necessary to know | 
the slope of the free energy curve i.e. Of/on. This is most conveniently 
expressed by considering n=«n, at any point on the curve. Thus eqn. (7) 
becomes 


of aN 
iP el ay ag 
N an —KT' In ome +W+kT In C—kT In 2. 
Now N-—on=N—n, since n<N and ny <N. 
Hence 
Clara x a6 2 : 
N an =khT In «+ E In yar +WHkT In C—kT In 2 | : 


But the term in brackets is zero, therefore 


0 
v2 ering. See eee En VEY 
3.3.1. Free Energy Curve for the Liquid Phase 

Unfortunately, the form of the free energy curve for the liquid phase 
is not known. Nor are the concentrations of copper in the liquid phase 
small. Thus the liquid cannot be considered a dilute mixture. However, 
it will be seen that if it be assumed that the liquid is ideal then the 
subsequent calculations agree fairly well with experiment. 

For an ideal mixture, the heat of mixing will be zero, and the entropy 
of mixing ideal (Lumsden 1952). 
Let 

F (,°=free energy per atom of pure liquid copper at 7'°K. 


F',,°=free energy per atom of pure liquid germanium at 7'°K. 


Then the free energy per atom of a mixture of (V—n) germanium 
atoms with n copper atoms will be 


l kT 
F= V [(N—n)F go+nF oy?) + V (N—n) In (N—n)+nIn n—N In v | 
OF , = n 
therefore N a =(Loy°—F go) +kT In EES 


=3P+i In(s"-). ia tee 
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3.4. Practical Application 

3.4.1. Physical Picture 
__ Before proceeding to the numerical calculation of the solid solubility 

curve of copper in germanium, a simple picture can be given which 
explains the shape of the curve, and gives physical meaning to the 
calculations. The solid solution will be referred to as the S-phase. 

Figure 3 shows schematically the free energy curves of the liquid and. 

f-phases at a fixed temperature 7 just above 650°c. Curve 1 is the free 
energy curve of the liquid phase, curve 2 that of the B-phase. The right- 
hand side of the diagram represents 100% Ge. 


Fig. 3 
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—~< COPPER CONCENTRATION 100 °/. Ge 


Variation of free energy of 6 and liquid phases with composition 
(T just >650°c). 

The free energy curve of the B-phase is shown rising sharply as the 
concentration of copper tends to zero. This is because the configurational 
entropy term in the free energy decreases rapidly in this region. 

From the principle of lowest free energy, the boundary points of the 
liquid and f-phases are given by the compositions at A and B respectively. 
In fig. 3, the tangent AB is drawn nearly horizontally, and the point B 
is thus nearly at the minimum of the free energy curve of the f-phase. 
If it is assumed that below a certain temperature this condition holds, 
then it would be expected that in this region 2,==% (see §3.3) and that 


W 
n.=-4 exp (— zp) oe, ee, 


where n, is the saturation concentration of copper atoms in the f-phase, 
i.e. the same form of temperature dependence is obtained as was found 
for thermal acceptors between 550°c and 800°C 
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As the temperature is raised curve | will sink relative to curve 2, and 
there will be a tendency for the point B to move to the right of 
the minimum. This tendency will be small until a limiting temperature 
is reached. Above this temperature, n, will deviate strongly from the 
value given by eqn. (11), and eventually will decrease with rising 
temperature. This state of affairs is shown in fig. 4, where the notation 
is the same as in fig. 3. As the melting point of pure germanium is 
reached, the two points A and B will coincide at 100% Ge, i.e. the 
solubility will have become zero. 


Fig. 4 
‘curve 2 
‘Ny () 
a 
CURVE | 
(LIQUID) 
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FREE 
(encrev 


LIQUID 


COPPER CONCENTRATION 1O0%Ge 
Variation of free energy of 8 and liquid phases with composition (high temps.). 


3.4.2. Numerical Calculation 
Now for a given temperature 7’, 


sh _ (4 
dn n=n, dn N= Ng 


(see § 3.2) i.e. using eqns. (9) and (10) 
kT In a=6F9+kT In he 


N—n, 
Nien 0 
Hence ‘ a = : exp Ga : 
But 2 =o 
therefore (- =") Ng—=Ng EXP Ga 


2 N=n\ sox, (SPO—W 
.e. My N,=A exp ie. ~ ae . (12) 
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By considering the quantities in this equation in terms of partial thermo- 
dynamic quantities, it can be shown to be equivalent to eqn. (2) apart 
from the term y,! since 


xEXPTL. Hs NTS 


3 . 
Plot of In (=) Ns against ~ (values gained from expt. points for nq). 
NL ; 


Now, 5F° would be expected to vary little with temperature, i.e. if 
In[{(N—n,)/n,}n,] is plotted against 1/7,a straight line of slope (8f°— W)/k 
should be obtained. These quantities have been plotted in fig. 5, where 
it can be seen that the linear relationship is obeyed fairly well, with 

6 F9— W=—2-12 ev 


and A= 1-91x 106/em*. 
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In fig. 6, the Inn,—1/T7' curve corresponding to this line has been 
plotted, and compared with the experimental points. 


Fig. 6 
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~ 
Comparison of: theoretical curve for thermal acceptor concentration with 
experimental values. 


§4. Discussion 


It must be emphasized that the experimental values of n, are con- 
centrations of thermal acceptors, and not explicitly of copper atoms. 
Recently, van der Maesen and Brenkman (1954) have shown that nickel 
also diffuses rapidly in germanium. Their measurements of its solid 


solubility indicate that it is smaller than that of copper by a factor of 


10 at 700°C, changing to a factor of 4 at 920°c. The effect of the presence 
of nickel on the measured values of n, would thus be small except at the 
higher temperatures. The tendency of the plot of In [{(V¥—n,)/n;}n,| 
against 1/7’ to curve upwards at higher temperatures (fig. 5) may possibly 
be explained by the presence of nickel. However, there is good agree- 
ment between the values of n, used in this paper, and the concentration 
of thermal acceptors found by Fuller and Struthers (1952). These, in 
turn, were in good agreement with the concentration of copper atoms 
diffused into the identical samples. 

The variation of the solid solubility of nickel (van der Maesen and 
Brenkman 1954) with temperature follows the predictions made here, 
although the exponential low temperature region is small (above the 


eutectic temperature). The measurements also agree well with the type of 


variation given by eqn. (12). 


——— oe 
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LI. A New Treatment of Anharmonicity in Lattice Thermodynamics : I 


By D. J. Hooton 
Department of Mathematical Physics, University of Edinburgh * 
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SUMMARY 


Born has given a method by which the anharmonic vibrational motion of 
the atoms in a crystal can be approximated in terms of an adapted set of 
harmonic oscillations, these differing from the usual modes of vibration of 
harmonic lattice dynamics ; this method is here redeveloped from another 
standpoint and extended to give explicit results. The thermodynamical 
formulae of the anharmonic crystal can be given a simple form: they 
consist of the customary formulae for a set of harmonic oscillators, but 
now expressed in terms of new frequencies defined essentially from the 
quadratic and quartic terms in an expansion of the potential energy, plus 
correction terms which express the difference between the actual anhar- 
monic motion and the harmonic approximation. For small anharmonicity 
the formulae reduce to those of the usual perturbation procedure (allow- 
ance for thermal expansion) with some extensions, but they also give a 
solution in the case of strong anharmonicity. This latter solution will 
later be used in a discussion of solid helium. 


§ 1. 
In an attempt to establish the existence of intrinsic imperfections in a 
crystal at thermal equilibrium on the basis of lattice dynamics Born 
(1951 a) has paid critical attention to the anharmonicity of the atomic 
vibrations in a way which allows quite another type of problem to be dealt 
with. When the vibrational energy is large compared with the static 
cohesive energy the displacements of the atoms, say from regular lattice 
sites, may far exceed the limits of smallness implied by the use of the usual 
harmonic theory. For example, Domb (1952) has shown that in solid 
helium the zero-point energy (some 10-100 times larger than the thermal 
energy in the region 0—20°R) is such as to produce relative displacements of 
an order nearly three times as great as normally found (for example, in 
the other inert gas solids). Corresponding to this large energy, essentially 
anharmonic in nature, and one and a half times as great as the lattice 
potential, is a large vibrational expansion of the lattice as a whole. The 
usual harmonic lattice dynamics is characterized by a frequency propor- 
tional to »/(d?U/dV?), where U(V) is the static lattice potential, say per 


* Communicated by Professor M. Born, F.R.S. 
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mole, and V the molar volume; it is found for helium that d2U /dV* goes 
over to negative values well within the experimental range of V so that 
the quadratic theory is inapplicable. Clearly, to describe such a solid we 
- must pay proper attention to the anharmonic terms in the expansion of U 
and Born’s considerations provide a direct means of approach. A similar 
situation may be met with near melting in an ordinary substance. 

In the present paper general formulae will be developed by an alternative 
method to that of Born; the assumptions inherent in his treatment can 
then be clearly seen and it is possible to extend his results and to obtain 
explicit solutions. By suitable approximations these can be worked out 
to give a reasonable description of the properties of solid helium, and this 
application will be made in a later paper. 

In general, the potential energy U of a vibrating molecular system is to 
be expressed as a polynomial expansion U(x) in the individual atomic 
displacements x,—2x,°, with suitable reference positions z,°. The first 
point of Born’s argument is to avoid the usual static condition dU/dx,°=0 
as a definition* of the x,°, since the mean atomic positions (x,) diverge 
increasingly from these static x,° with increasing distance through the 
solid from any atom chosen as a fixed origin. Instead, the x,° are defined 
as the <(a,,) themselves; this preserves the idea of relatively small dis- 
placements x,—7x,° so that the expansion of U(«) is justified. The new 
definition may be expressed in terms of the first and third degree terms U,, 
U, of the expansion taken together. 

Naturally, if the z,° are assumed to form an infinite lattice structure at 
the beginning the individual static conditions 0U/dx,°=0 define only the 
relative basis vectors, being otherwise identically satisfied ; they are in 
effect replaced by conditions 0U/da,,—=0 derived from the static part of the 
free energy (free crystal; a,,—a, .a,, the a, being the cell vectors, «, B 
=1, 2,3). Born’s expansion difficulty is avoided by the usual procedure 
of including the vibrational terms of the free energy (~U 2) in the definition 
(thermal expansion) ; the result is of course exactly equivalent to that of 
the new argument, which is in this case superfluous. Its relevance is toa 
finite lattice or a general molecular structure, and since the equations are 
non-linear there is the possibility of some irregular lattice solution ; 
however, it now seems unlikely that such a solution, more stable than the 
regular lattice, will in fact exist. 

In the following we shall retain the individual coordinates 2 eeven 
though our interest is to develop a general description of strongly anhar- 
monic vibrations in a regular lattice ; the work is then symmetrical and no 
more complex, and may also be applied to other cases. For example, it is 
of some interest to study the linear monatomic chain, which can be solved 
explicitly ; the theory shows exactly the trivial nature of the modifications 
ment Mee ee 

* This definition is a necessary consequence of Born and Oppenheimer'’s 
general treatment (1927) of coupled atomic vibrations ; however, Born oe b) 
has given an elegant generalization in which the static condition is no longer 


necessary. 
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in a finite but otherwise regular structure which occur at the free surface— 
a particular justification for the cyclic boundary condition—and will be 
given elsewhere. 

The second, and in practice more important, aspect of Born’s method 
is based on the necessity (for example, in solid helium) of considering that 
the atomic vibrations may have amplitudes which, though relatively small, 
are not so small that the corresponding frequencies can be specified by the 
quadratic terms U, alone. The new definition of the 2° of course intro- 
duces a first change in the frequencies, arising from the third degree 
terms U;, but does not affect their functional form. Obviously fourth 
degree terms must be included to make any functional difference, and it is 
generally true that essentially the terms of odd degree in U(x) have to do 
with the structure of the crystal, the terms of even degree with its energy. 
The usual harmonic theory considers first and second degree terms to- 
gether; any extension should clearly include third and fourth degree 
terms together (as we shall do here), and so on. 

If the third and fourth degree terms in U(x) represent only small 
corrections to the linear and quadratic, the explicit form of Born’s 
results developed by the present author can be immediately expanded in a 
perturbation procedure (see Born 1952) ; this will be illustrated elsewhere 
for the linear chain. We may note here that the zero order solutions are 
simply the solutions of the usual harmonic theory, the equations 


dU jor P=0; 
va 0U7/dx,° ox, . Cry Cy p=; 85475 . . . . . (1) 
kk’ 


defining the x,° together with the frequencies w; and modes e,; of the 
proper vibrations. First order corrections 62,° arise directly from the 
terms 0?U/dx,° 0x,° 0x," of Us, while the simultaneous corrections 
dw;, 5€,; depend directly on the terms 04U/0x,° 0x, Oxy OX, of Uy 
and are not given solely as the result of thermal (anharmonic) expansion 
expressed through the dx,°. 

The nature of Born’s approximation can be illustrated by the following 
simple example. In the problem of the vibrations of a molecule, one has to 
solve Schrédinger’s equation for a potential well of shape (a). 


—_ 
=e 


Shape (a). Shape (0). 
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When the vibrations are small (energy level as shown) only the bottom 
of the well comes in and the harmonic approximation (the dotted curve) 
is reasonable. If, however, the well is of a different shape (shape (b)) 
(consider, for example, an atom in the interior of a crystal, with large 
vibrational energy) this approximation breaks down completely (dotted 
curve). Nevertheless, a pseudo-harmonic approximation is still possible 
with a harmonic oscillator wave-function and the frequency parameter 
fitted by a variational principle equivalent to the fixing of average position 
and mean-square deviation mentioned above (dashed curve). Since the 
energy is always large in this case the central hump is of no consequence for 
the pseudo-harmonic vibration. 


§ 2 
Ss <- 
Following Born we write the Hamiltonian of a crystal containing 
NN atoms (k=1...3N) in the form 
Mg rie EA sO (Ls, Kan eeee en (ae ee es 12) 
i 


with 
U(x)=U,+U,(x)+U,(4)+ ..., 

1 
U,(x)=— 2 Une: ME LN Le) ny OU ey ee =O Ul Or OG, Oo (3) 

“(k 
the sum being over all sets of r indices k, and the 2,° having arbitrary 
values. As in the usual theory of normal coordinates, we make an 
orthogonal transformation of the form 

| 


1 

: = x 

Uy — XL yp = TF = C75 p> 

VM, j 
Po A Mh, XC Djs (4) 

a 

2 Ong Cp =9 597s 21 5 Cnjg=O nk’ 


but now with arbitrary values for the e,;. The potential U(x) then goes 
over into a polynomial expansion 


0(q)= Uo +01(G)+%2(9)+ -- +; 7 
1 hos Ong Cnij? ++ Tie 
v,(q)= al oor “arts 9595" ELS Viy" ‘ — Jy My ss kh’ | 
(the sums being over all sets of r indices k, j) ; the energy (2) becomes 
H=t J pP+v(q). aa (Wires tse fesse (20) 
j 


In the usual theory the normal coordinates q; are chosen by making the 
quadratic term v,(q) diagonal, and the motion 1s supposed to be one of 
small vibrations about the positions of static equilibrium : 


v(Q=t 2 0? G7; v1(q)=9. eo +(6) 
j 


The first condition defines the modes of vibration C xj together pa 
harmonic frequencies w,, exactly as in (1), while the second defines the x, 
according to the static conditions U ,=0, also contained in (1). Generally 


195 2G. 
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the higher degree terms of v(q) are ignored in (2’) or they may be treated as 
a small perturbation once the modes of vibration have been defined. 

In the present theory we retain the transformation (4) but drop the 
subsequent conditions (6); further, the anharmonic terms of v(q) (at least 
up to v,(q)) are considered an essential part of the Hamiltonian of the 
actual crystal. This Hamiltonian (2’) can now be written out as 


=% x pP+ Uo+01(q)+v2(¢)+s(q)+e(g)+---> + + - (7) 


where v,(q) is no longer diagonal. Born’s method is to compare this 
energy with that of an arbitrary set of harmonic oscillators, 


H=} 2X pP+Uot+3 & w;?? g;?, 1s WORE AG 
) j 


new frequencies w,, together with new transformation parameters e;,; and 
new vibrational centres 2°, being chosen to make (8) the best possible 
approximation to (7). From (7), (8) 
H=H-+H’, | 
H’=v,(q)+(v.(q)—4 a w; q;") +03(9) + a(Q) + ee J 
j 


the statistical problem involved is to find the modification in, say, the free 
energy F (which corresponds to H) due to the perturbation H’ under a 
suitable choice of the parameters w,, ¢;,;, v,°. This choice must in the 
first place allow such a perturbation to be made, so that the changes L,’, 
say, in the energy levels #,, of H are small by comparison with the £Z,, 
themselves. It must also agree with the condition <#)=x° andy if the 
anharmonicity is neglected, lead to the customary w,;, €,;, #,,° defined in 
(1) or (6). 

Instead of considering with Born the statistical density matrix p(q, q’) 
we shall work from its integral, the partition function 


(9) 


- ) 
F=| p(q, 7) dq=~ exp (—BE,), 
|e ee aaa! (10) 
6 ,=L,+£,, B= VkT. 
Script letters will be used throughout to refer to the actual Hamiltonian 
(7). The partition function of the system of oscillators, together with the 
thermal average of the H,,’, taken over the unperturbed energy levels £Z,,, 
will be written : 
Z=Lexp(—BH,), H'=Z-12'H,' exp(—BE,);. . (10’) 
nN n 


otherwise a bar will be used to denote the average value of any function of 
n, formed in respect to the partition function Z. Thus (10) becomes 


%=Z exp (—BE,’)=Z exp (—BE’) exp[—B(E,/—E)]. (14) 
Supposing an expression can be found for the final term in (11), say* 
exp [—f (£,,—E’)|=exp(—fI), . . ae (12) 


is More particularly, we could write etl, since the odd terms in the 
expansion of (12) disappear on averaging. 
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we have at once for the free energy ¥ of the anharmonic crystal 


A=—kl log S=— kl ioe ZH AT (13): 
- where it is well known that for the Hamiltonian (8) 
—kT log Z=U,+kT 2 log [2 sinh Gohl): vase oy (145 
j 


The thermodynamical properties of the anharmonic crystal can therefore 
not be described by an ‘adapted’ harmonic theory (with frequencies 
fixed anharmonically) alone; there remain the additional terms in (13), 
which represent a correlation between the auxiliary variables q;- These 
terms, as well as the harmonic part of ¥, are of course functions of the 
parameters w,;, €,;, 7,° and of the temperature. 

Since for a lattice the «,° replace the usual macroscopic parameters @,, 
defined in §1, the last remark suggests immediately that they should be 
chosen (simultaneously with the w, and ¢,.,) to give a minimum value of F, 
in analogy to the minimum conditions usually imposed macroscopically. 
It is easily shown that to the 3N coordinates x,°, the 3N frequencies w,; 
and the $3N(3N —1) independent quantities e,, there correspond just so: 
many independent minimum conditions 


OF 7 

— = k—1—3N, | 

5,0 = 9 | SN, | 

OF : 

—o == (), j=1—3N, | * A ° (15) 
OF OF 


= ——- Rt nA | 
7 ee ; ens) 220.4555) ,—=1-- 3 


k je J 
by which to determine them. If there are external stresses they must of 
course be allowed for by writing the equivalent forces —f, on the right 
hand side of the first set of equations in (15). For a regular lattice these 
equations are equivalent (apart from a basis) to 0.7 /da,,=0, subject, say, 
to a possible condition dF /0dV=—p. The relation between (15) and the 
conditions introduced by Born is discussed below. 

The success of this formal development depends on our being able to 
calculate the average perturbation H’ and the quantity I’. It is clear that 
the energies &,, E,,, H,/ are all proportional to N, the size of the system, 
wereas 8-1—kT' corresponds to the energy of one degree of freedom alone ; 
thus to have any significance the exponentials in (11) should be written 


where 7,=H,'/N, y=H'/N are energy changes per oscillator. Thus for J” 
the problem is to find the averaged Nth power of the exponential 
exp Le (7,—n)|=exp [—B (7, —n" )J=1—B (7) n 7 )+ sinees)|9) <* (17) 


where we may split off the change in the zero energy Ky so that 
7p =(E, —Hy )/N. The result itself should be put in the-form of an 
Nth power since its logarithm appears m the free energy (13), and the 


2G2 
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latter is itself of order VN. We have essentially to do with an expansion in 
powers of 8, and the convergence depends on the order of magnitude of 
the thermal energy. For example, in solid helium the vibrational energy 
is ~10° erg/mol., the thermal energy 107-10*% erg/mol.; at the solid 
helium temperatures in which we shall later be interested (~10°K) 6 is 
~10% erg-1, so that 7,7B~10-! (N~10?3/mol.). Thus at such low 
temperatures an expansion in powers of 8 converges slowly and we can not 
expect that an accurate evaluation of I’ will be easily made. Instead we 
shall make a gross approximation by taking only the constant term in (17), 
which implies that [=0: even for helium, however, this gives reasonable 
results. 

To evaluate the remaining quantity #’ we must express the various #,,’ 
by the usual perturbation formulae in terms of the matrix elements 
H,,,,, of H’ with respect to the eigenstates of the oscillators (8). Going 
back to (9) we see that H’ consists of a sum of homogeneous polynomial 


/ 


expressions in the oscillator amplitudes g;. The H,,,,,/ may be written 


down at once from the well-known values of (9,4; - « -),,., and from them 
the energy perturbations H,’, H,’? ,..., and so the thermal averages 
H'), H'?) ,..., in successive orders. Thus, 
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n 
all odd order terms vanish by symmetry, and we have used the fact that 
qji=39,;?2. For the (non-degenerate) zero-energy perturbation Hy’ (which 
is in fact all that is needed in the helium problem to be discussed later) 
it is found that in the second order 
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7-0 
But since #y= lim #’® we may infer that H’® is to be got from (19) 


an: T-0 . . . 
by writing q,? instead of its limit a,?.__The general form to be expected in 
the higher orders is then clear from (18), (19): appropriate powers of the 
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curly brackets plus some remainder terms. The latter diminish rapidly 
in magnitude (see below) and as a further simplification we shall neglect 
them in the second and higher orders. 

If now (14), (18) and the generalized form of (19) are put into (13), 
together with the approximation [=0, and the minimum conditions (15) 
applied, some lengthy but straightforward calculation shows that a 
solution can be obtained by setting the curly brackets equal to zero : 


: = (20) 
Vip FEZ Vig p50" Gr=0 5 


here terms of the fifth degree in ¢ are dropped, since we are considering 
only a first correction (terms v3, v4) in the anharmonicity. The eqns. (20) 
for the parameters x,°, w;, €,; also ensure that we have a proper pertur- 
bation procedure : the brackets in (18), (19) and the higher orders vanish, 
and since the average effect of fourth degree terms is in general to be 
treated as equivalent to that of the quadratic, the order of magnitude of 
the remainder terms, for example in Hp, is given by 
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and similar fractional expressions. Transforming (20) according to the 
reciprocal of (4) we get 
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an immediate extension of the old definitions (1). Inthe case when U,, U, 
are small compared to U,, U, (21) may obviously be solved by successive 


approximations. 
Just these relations (21) have also been found by the present author as 


the explicit expression of the general conditions originally developed by 
Born (1951 a) from the equations 


{qj )=9, or {U4 )=X;,", (G4 =4%2 Oy 415 ase ARF 


as definitions of the parameters x,°, w;, @,;3 « ) denotes the thermal 
average under the Hamiltonian (7). Born’s development is based on an 
expansion of the density matrix p in (10) according to the expression 


p=99, g=14+-9179et --ss 
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where Z= fo(q, q) dq, and g, is a polynomial of degree r in the coordinates 
q; by comparison with (11) it is clear that 


1+9,+9ot ...=exp (—BE’) exp (—BP)=(1—BE"— .. .} 


(1446? (H,’ B+ ...). 


Thus the polynomials g,, like the #’, are of the order V, and to make a 
properly convergent expansion we should have to write 


p=e(g')", gl orgs ae ee 


which destroys the relative simplicity of Born’s equations. However, 
as his method stands, the expression for g got from the first few polynomials 
g, turns out to be the same as the leading terms of 1—BH’™. No higher 
powers of f are obtained before the work becomes prohibitive, so that in 
practice the term in I’, which requires the use of (23), is automatically 
excluded, nor is there any contribution from the #’@ or higher orders. 
Born’s treatment by successive polynomials of course implies that the 
vq) are of successive orders of magnitude and so gives the terms of 
1—fE’™ only in successive approximations. The present treatment is 
not so restricted and allows, in particular, v, to be equivalent in its effect 
tov,; the simplifications which exclude the higher orders and the term in J* 
are thus of more general significance than Born’s method would imply. 
The approximate identity between the results of (15) and (22) confirms 
that the present method of specifying the parameters 2°, w;. ¢,; has a 
good physical meaning as discussed in § 1 above. 


§ 3. 

A solution of (21) in which the anharmonic terms are not regarded as 
small corrections is evidently connected with the presence of large vibra- 
tional energies; the determination of the frequencies is thus the main 
problem and the exact form of the crystal (7,°) and of the modes of 
vibration (¢,;) is relatively unimportant. Whatever the «,°, we may 
isolate the effect of the new vibrational energies (frequencies) from the 
second set of equations in (21) by assuming that the e xj have the customary 
form defined by the non-diagonal equations (j4j’)in(1). Of course this 
implies that Uj, has been made diagonal, and we now write the diagonal 
equations of (1) as 


yw e€ kj e€ kj 
eed | 
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to avoid confusion with the new frequencies w;. The assumption refers 
only to the non-diagonal equations in (21): 
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and the approximation implies that the last term in the bracket (though 
not small) is largely independent of k, k’. The diagonal equations remain, 
and may now be written 
wP=w +t 2 Vy 5379 Ws Ge lIeBN. se oe. (25) 
J 
Although there is no necessity for the w, to be positive (roughly speaking, 
they are positive only for atomic spacings less than the inflexion distance 
of the interatomic potential) that no longer prevents a description of the 
vibrational energy being made, so long as the whole expression (25) leads 
to positive w,?. 
The eqns. (25) are non-linear, since 


G2=E,|w?, | (26) 

B,=Yie, coth (HofkT), [f° 9 °° °° © 
the latter being the thermal average of the energy of the jth oscillator 
alone. Thus (25) must be solved by iteration or approximation. At low 
temperatures, or for large zero-point energy, the use of a,” as in (19) leads to 


ww, + dh 2 Vyyq79°|@y > . . . . é ; (27) 
4) 


the equations in this form can be solved explicitly for the linear chain, 
and the resulting solutions will later be used to give an approximation by 
which the helium problem can be discussed. 
Going back to (13) we see that, subject to the approximations that have 
been made, the free energy of the anharmonic crystal is to be written 
F=U,+F-4 Ff, E, | 
#) 
fy=1—0;” w;”; | 
F=2,F, now refers to (14), or (8), without the constant (7), while the final 
term of F expresses the remainder of (18) with due regard to (20), (25). 
Some care is needed in general in applying the formulae of thermo- 
dynamics to (28) since the new frequencies w,, solutions of (25), are expli- 
citly 7'-dependent as well as being functions of the x,°; the w,° and ¢,; 
are functions of the x,° alone by (1) or (1’). In a lattice the x,° are 
equivalent to the macroscopic lattice parameters, and in any case they 
specify the molar volume V ; thus, excluding any relative variation of the 
2,2 with T at constant V, we may write 
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For example, the entropy, energy, and specific heat at constant volume 
follow from (28) according to 
S=S+1ELX Ff, Cpt 2'8;, 
j j 
é= U,+H#+t xf; (CL Ea a ie ee ee (29) 
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the final terms contain all the variation of w at constant V, while the 
S=2,8,, H=2,E;, C=2,C, are the usual oscillator functions, It is to be 
noted that the latter are derived from (28) by the operator (0/07’),.7 
alone: operations at constant volume in the ordinary oscillator thermo- 
dynamics imply that the frequencies are held constant, and corresponding 
oscillator relations must be used here in the same sense. 

In practice it is to be expected that the anharmonic properties of the 
crystal will depend mainly on volume, so that the temperature dependence 
of the energy (frequencies) at constant volume is a small correction ; 
that is certainly the case for solid helium, for which the eqn. (27) may be 
used, and the terms in s,, e;, c; may then be safely ignored. 

The whole problem is now to describe the new vibrational spectrum w,, 
and from it the f; of (28), in any particular example. By (25) or (27) this 
spectrum is to be got from the fourth degree terms in the expansion of 
U(x) as well as from the quadratic and the solution can not be simple. 
For practical purposes we might hope to make an adaptation of the well 
known Debye approximation and this will be discussed in a following paper. 


ACKNOWLEDGMENTS 


I should like to thank Professor Max Born for suggesting this problem 
to me and for continued helpful discussion throughout. The work was 
supported by a Mr. and Mrs. John Jaffe Donation from the Royal Society. 
I am also indebted to the Referees for helpful criticism and comment, and 
in particular for the formulation of the example given in the first section. 


REFERENCES 


Born, M., 195la, Fest. d. Akad. Wiss. Gottingen (Math.-Phys. K1.), ‘‘ Die 
Giiltigkeitsgrenze der Theorie der idealen Kristalle und ihre Uber- 
windung”’; 1951 b, Nach. d. Akad. Wiss. Géttingen (Math.-Phys. K1.), 
“ Kopplung der Elektronen- und Kernbewegung in Molekeln und Krist- 
allen’; 1952, Comptes Rendus 2° Réunion Chimie Physique. Paris 

Born, M., and OppENHEIMER, R., 1927, Ann. Phys., 84, 457. 

Doms, C., 1952, Comptes Rendus 2° Réunion Chimie Physique. Paris. 


c eR 


LU. A New Treatment of Anharmonicity in Lattice Thermodynamics : IT 


By D. J. Hooton 
Department of Mathematical Physics, University of Edinburgh* 


[Received November 6, 1954] 


SUMMARY 


The thermodynamical formulae of the previous paper are worked 
out with the help of an adaptation of Debye’s continuum approximation ; 
in particular, the specific heat at constant volume is put into a form 
suitable for numerical calculation. This formula contains, however, a 
factor which expresses the (possibly strong) volume dependence of the 
relation between the new frequency spectrum and that of the customary 
lattice dynamics ; the factor appears in addition to the Debye character- 
istic temperature 6 and must be estimated in any particular application— 
for example, in the following consideration of solid helium it will be 
approximated from a linear chain model. The meaning of a Debye 
characteristic temperature in the anharmonic theory is discussed, and 
the place of an empirical Debye temperature, determined by fitting 
specific heat measurements to a theoretical specific heat formula, is also 
considered. A discussion of this fitted Debye temperature (due to 
Domb and Salter) is adapted to the anharmonic theory in order to give 
later a correct application to solid helium. 


$1 

In a previous paper (referred to as I) a method originally due to Born 
for treating anharmonicity in a crystal has been given explicit form in a 
way suited to the discussion of problems involving large vibrational 
energy. In this paper we shall deal with the application of Debye’s 
continuum approximation in the new theory. The theoretical definition 
of the Debye characteristic temperature @ gives a result of the same form 
as proposed by Domb (1952) in connection with solid helium by fitting 
the Debye specific heat formula to experimental measurements. Such 
an empirical parameter, whose variation at constant volume is useful 
as a guide to the exact lattice spectrum, can also be defined in the new 
theory, and the justification of Domb’s formula will be discussed. 

The results obtained here will be illustrated elsewhere for the example of 
a linear monatomic chain, which can be solved. explicitly ; they also 
provide a basis for the thermodynamical discussion of solid helium to be 


made in a following paper. 


* Communicated by Professor M. Born, F.B.8. 
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§2 

The w,° of I are assumed to form a regular cyclic lattice in what follows, 
and for simplicity we shall deal with a monatomic lattice only; there 
would be no essential difficulty in treating ‘ optical ’ vibrations in addition _ 
to the ‘ acoustic’ vibrations discussed here. With the 2,° replaced by 
lattice parameters a,,, the first set of equations (I, eqn. (21)) is to be 
replaced by conditions* 0F/da,,—0 (free crystal; otherwise external 
‘forces’ must be written on the right-hand side, for example—and this 
is important for helium—d¥/0V=—p). The remaining equations 
(I, eqn. (21)), approximated according to (I, eqn. (24)), namely 
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give simultaneously new harmonic frequencies w; which are functions 


of the a,, and of 7’, while the w,® and e,; are the well-known harmonic 
lattice frequencies and plane-wave vectors, functions of the a,, alone. 
By (1) there is a 1: | correspondence between the w; and the w,° and it 
may be shown that w; and w;° vanish together; we may therefore speak 
of three acoustic branches to the new frequency spectrum. However, 
although the corresponding vibrations have the usual form (in this 
approximation), there are no longer any linear wave equations by which 
to determine their frequencies. 

Failing an exact solution of the non-linear eqns. (1) we can apply the 
well-known Debye approximation, most simply in the isotropic form. 


Thus the density of vibrations in the continuous Debye spectrum, 


9N 9NT 1 
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replaces the unknown dz of (1). Here the cut-off frequency w”, or the 
characteristic temperature 6, is strongly volume-dependent, with in 
general a 7'-dependence at constant volume ;{ the latter is to be ignored 
when the corresponding property of the w,; is ignored, as mentioned at 
the end of I. 

The thermodynamical functions (I, eqns. (28), (29)) of the anharmonic 
lattice depend on the frequencies w, through the usual oscillator functions 
derived from I’, (the free energy), but as well contain additional terms 


* See, for example, M. Born (1939). 

+ This of course has nothing to do with the empirical 7'-dependence at 
constant volume of a * fitted’ 6, got from specific heat values, which expresses 
the ‘ non-Debye ’ nature of the actual spectrum. 
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which provide corrections to the new harmonic approximation. Thus, 
for example, 


E€=U y+ E—45f,(Ej;—TC,)+5¢,, 
; 
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Ow; ) : : : r io) 
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where H; is the energy and C;, the specific heat at constant volume of an 
oscillator of frequency w,, as defined in the usual theory. 

The Debye approximation (2) can only be applied directly to functions 
of the frequencies w,; which do not depend otherwise on the index j, and 
this excludes terms containing f; or the derivatives of w, such as appear 
in (3). It is not difficult to show that all the terms in (I, eqns. (28), (29)) 
which are not functions of the w,; alone are of the form 27; 42;, where 
92, 1s such a function, and 7; is equal to (dw,/0T')y, (dw,/0T)p?, (0?w,/0T?) p. 
the product of these with f;, or f, itself. To get any further we must 
replace these 7; by quantities which are constant along the spectrum, 
which is similar to Griineisen’s approximation, that all dw,°/dV are the 
same, in treating thermal expansion. Strictly speaking, we should 
replace the vibrational index 7 by the reciprocal lattice indices 


(), h=1—N, j=1-—3 


(acoustic branches); the 7; might then be averaged for each branch of 
the spectrum separately. However, to use the isotropic Debye density 
(2) we need overall values of the 7 and need not introduce the separate 
indices. So far as f, is concerned, the assumption that w,;?/w,? is roughly 
constant along each branch is not an unlikely one, it being shown else- 
where that in the linear chain model f,; 7s constant along the spectrum. 
We have already noted that the 7-dependence of the w, is most probably 
a small effect, if not to be ignored, so the averaging of the derivatives 
(dw,/6T')y is unimportant. 
ON 713 “O/T 
Thus, 27; 82; = ee 7 | Q(E)\2 dé, 
j 0 

in which 7 is strongly volume-dependent and has also, in general, a small 
T-dependence at constant volume. The formulae (I, eqns. (28), (29)) 
then go over into Debye integrals such as (from (3)) 
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which are readily integrated in the usual way. For the application to: 
solid helium we shall drop the terms in (dw/07')y throughout. In that 
case (4) becomes 


E=Uyt(1—df)EP(6/T)+4fTOM(H/T); . . - - (8) 


we may also write down the specific heat at constant volume from 
(I, eqn. (29)), or by differentiating (5) : 


= (1+4/)0(G/7)— 2/08 (6/1). eee 


Here HE”, C” are the usual Debye formulae, and C” the Einstein formula 
for the molar specific heat, defined as follows : 


B(6/T)=3RTD(6/T),  C(6/T)=3R{4D(6/T)—3P(6/T)}, | 
C*(6/T)=3RQ(6/T) ; 


| 
[3 polT 
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We see from (6) that in this adaptation of the Debye theory there is no: 
longer a universal specific heat function in the variable 6/7’ alone; an. 
explicit volume dependence remains through the quantity /. 

The problem of finding the new harmonic spectrum in any particular 
example is now replaced by that of determining @ and the quantities 
Peco OL! ersten Of the latter, only a suitable average f is of real 
importance, and it may be possible to simplify its determination by the 
use of symmetry properties of the lattice—but now the symmetry of 
fourth degree terms is in question rather than that of the usual quadratic 
terms alone. In applying the theory to solid helium we shall use as an 
approximation to f the value derived for a linear chain of helium atoms ; 
this exhibits the strong volume dependence which is the new factor in: 
the theory. 

For 6 we might follow the empirical practice of fitting the specific heat 
formula of the Debye approximation—now given essentially by (6)—to. 
experimental values ; it is necessary to distinguish the resulting ‘ fitted ’ 
parameter @ from what would be got by the uswal empirical fitting of the: 
Debye formula C? alone. The latter, say 6’, is discussed for helium by 
Domb (1952) and by Dugdale and Simon (1953) but has no real place in. 
the anharmonic theory since the properties of a strongly vibrational 
solid can not be regarded as those of a set of harmonic oscillators alone. 
# and 6’, as empirical parameters, are related by the equation 


€ (O/T) = CP (GPL) | te on 
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§3 

From a theoretical standpoint Debye’s approximation must be made 
so that the long elastic waves which it introduces are identical with the 
- long (acoustic) lattice vibrations, clearly elastic in character, that they 
have to replace. For the present Debye theory the harmonic ‘ lattice 
vibrations ’ in question are fictitious, being those of Born’s adapted 
harmonic approximation alone; they are not in themselves the actual 
(anharmonic) lattice vibrations of the crystal, and the intrinsic energy 
of the Debye continuum is not the same as the internal energy & given 
by (3). The parameter @ is in fact to be defined through the velocity 
of the plane waves e,; of the new harmonic approximation in the long 
wave limit j7—> 0 (or the averaged velocity over the three branches in 


the limit (‘) = (‘) ile, 3)e 
J 4} 


h » 2h 3f/3N 
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(9) is the usual relation between 6 and the average isotropic elastic 
velocity c of the Debye continuum, and we write c, for the limiting 
‘lattice ’ velocity of the harmonic approximation. In the usual theory c 
can be made equal to c, by a comparison between the linear wave equations 
of the continuum and those of the lattice; that is not possible here, 
since we no longer have any linear wave equations in the frequencies w,. 
To find c, it is necessary to solve the w,-spectrum itself near the origin 
of the wave-number space, and an explicit example of this calculation is 
given elsewhere for the linear chain. 

Of course, when macroscopic (i.e. experimental) elastic properties of 
the crystal are defined in the usual way—for example, an adiabatic 
compressibility «—they must contain all effects of the actual state of 
anharmonic vibrational motion, whatever it may be. Nevertheless, their 
definition implies that there are at least long elastic-type waves actually 
present which are harmonic in character ; this is not unreasonable, since 
most of the vibrational energy lies in the short-wave region, outside which 
the anharmonicity can be neglected. The corresponding macroscopic 
elastic velocity is essentially 


V Or 
(eK) — ama (sr) : . . . . . (10) 


(F would replace & for isothermal compressions; (10) is an isotropic 
average.) But it is clear that this concept of an effective harmonic 
continuum can not be extended to the short-wave region, and so used 
to define a Debye approximation (which refers to all vibrations), since 
the microscopic, and more energetic, part of the lattice motion must 
be regarded as anharmonic; we can see that harmonic Debye formulae 
alone can not give a sufficient description of the whole anharmonic 
motion since harmonic lattice formulae are not of themselves sufficient, 
even when the frequencies are anharmonically defined. 
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The present Debye velocity therefore satisfies a relation of the form* 


(10), V 07 il 


(a relation of this form must be true in any harmonic continuum), but 
with the elastic energy « equated to the energy U,+H of the new harmonic 
approximation itself. This energy is the energy of a complete set of 
harmonic vibrations, and differs from & by the final terms of (3). In the 
anharmonic theory the vibrational part H of « is essential; the equi- 
valence of c and c, then implies that the energy of the long waves e,,; (j — 0) 
is strongly dependent on that of the short waves, and such a relation is in 
fact contained in the eqns. (1). In the usual harmonic theory (relatively 
small vibrational energy) only U,) appears; the elastic velocity (11) is 


then 
V 02U 
Co ~ in JGR). ade a 


in immediate agreement with the velocity c, of the usual plane waves: 
of frequencies w,°. If correction terms H(w,°) are included, the change 
dc, corresponds to the dw,” given by the equations (I, eqn. (21)). 

From (9), (11) the volume dependence of @ is clear; there is also a 
temperature dependence at constant volume, through the energy term Z, 
and this has already been mentioned as representing the corresponding 
property of the w,-spectrum. Compared to the variation of 6 with 
volume (for example, in solid helium @ ranges from 40-160°K when V 
goes from 20-10 cm3), the temperature dependence is unimportant ; 
whenever (1) is approximated by the zero-vibrational form (I, eqn. (27)) 
E is to be replaced by the zero energy #, in (11) and this (theoretical) 
temperature dependence drops out. 

§, or c, is only formally determined by (11) since # itself depends on 
the unknown w,; and must be approximated by the Debye formula 
£°(6/T) ; thus we have an approximate differential equation for @, non- 
linear and of the second order : it is given essentially by the zero-energy 


form 
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* For this form, and for (9), see for example M. Born (1923). (10) corresponds 
to elastic deformations from a state of vanishing initial stresses, as is usual in 
elasticity theory. When, for example, an isotropic pressure is necessarily 
present—as it is for solid helium—we may expect a modification in the macro- 
scopic elastic velocity of the type 


V(8E/dV?) > (0? 6 /OV2 + fn (p/V)). 


It is reasonable that such a modification should also be made in (11); the 
energy of the anharmonic lattice under external stresses is then approximated. 
with the help of harmonic oscillations in a Debye continuum under ecorres- 
ponding stresses. 
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Such a relation can not be avoided in the strongly anharmonic theory, 
whereas in the usual case (with or without the correction 6Cy) it does 
not appear. <A formula of the type (13) has been postulated by Domb 
_ (1952) to account for the Debye parameter, called 6’ in §2, got by fitting 
Debye’s formula CO? as in (7) to the measured specific heat of solid helium. 
We shall next see how the analysis of Domb and Salter (1952) on which 
this is based is to be modified in the anharmonic theory, and whether 
it is of any significance for solid helium. 


§ 4 

Domb and Salter have considered the almost constant values of the 
Debye parameter 6’ found in many lattices at relatively high temperatures 
and any given volume when the specific heat formula C” is fitted to 
experimental results, and have given an expansion in powers of 


t=he,, [kT 


for the correspondingly small 7'-dependence (w,,, is the maximum lattice 
frequency); the variable x is ~ 6’/T and will in general be smaller than 
unity. In the same paper it is shown that the Debye zero-energy 2R0’, 
written in terms of the constant limiting value of 6’, is a good approxi- 
mation to the actual zero-energy of the lattice, which is rather surprising, 
since the consideration is based on (relatively) high temperatures while 
the zero-energy characterizes the opposite limit of large values of «a. 
We shall see below that in both the usual and the new theory the functional 
form of the constant limiting fitted Debye parameter as given by Domb 
and Salter’s method is just that of the theoretical expressions (12) or 
(13), as must be the case if this 0 or 6’ is to express correctly other thermo- 
dynamical. properties, such as the zero-energy, at different volumes ; 
only the constant factor is new. Moreover, the simple Debye formula 
2 RO’ as quoted above will be found to retain a certain validity in the new 
theory. 

6’ is defined by equating the Debye formula C”(6’/7)) to the experimental 
specific heat values. Instead of using an integral expression such as (7) 
for the molar specific heat, Domb and Salter use an expansion formula* 
due to Thirring : 


Sa her »p\2" } 
C= BNE {1+ 2(—I B55 Ter) } | 
v Sate (14) 
1 
Lp 
Dees Fe aN SOR | 


C! is the specific heat corresponding to any spectrum of independent 

3 . YD : A 

oscillators of frequencies w; and maximum frequency w,,. C” is got by 

a a 

* The factor 2n—1 is missing in Domb and Salter’s paper, although their 

numerical values are correct. Their formula for 6’, by comparison with (15), 
is also evidently in error. The B, are the Bernoulli numbers, 4, 3, -.--- F 
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writing in ,?, w? corresponding to the Debye spectrum; the experi- 
mental specific heats are represented when the proper p,;, w,, of the 
actual lattice spectrum is used. By equating the two expansions so found, 
and developing 6’ in powers of fiw,,°/kT’, Domb and Salter obtain the 
formula 


h 5? 1 /ho,,°\2 25 ee 
=a) | as eles — > mri 
v=7/( 3 )qr za of) (1, pi a') [72 teal 


the indices ° are introduced to correspond to the lattice spectrum w,° 
of the harmonic theory. Since p,° S «/m, where « is the force constant 
of the lattice, and with « S d7¢/da?, U»=N¢, V=Na’, M=Nm, the 
constant term of (15) is clearly* 


, 1 0h N13 V2 eu 


(16) is in immediate agreement with (9) and (12). 
In the anharmonic case the specific heat given in (I, eqn. (29)) is no 
longer represented by (14), but by 


onl )(led 2n 
@=3Nk{1+2(—1"B,—— 1, (Fe Loan 


w,; and w,, being the frequencies of the new harmonic approximation. 
Since 9’ is still defined by Debye’s formula, or (14), Domb and Salter’s 


method will now give 


Foe Bu 1 /hw,,\? 2(f—1) 
pal Cs—19) ao) 
25 (f—2) 
Wie eA a, 


‘On the other hand, if we fit the proper Debye temperature 6 of the new 
theory by equating @(6/7') from (6) to the lattice expression (17) we get 


ho | (5p 1 (he \? 2 f—1) 25 


(18) and (19) are in agreement with (8). The constant term of (19) is 
the direct extension of (16), but no longer refers to the empirical 6’ at 
all; the slightly different temperature dependences in (18), (19) are 
unimportant to the present discussion. Although in practice @ and 6’ 
are not very different, it is clear from the earlier discussion that the 
simple Debye formula C?(6’/7) has no real place in the anharmonic 


* Tf there is an external pressure p it is readily shown that 
Oh 7 U, ae 
da?” OV? 3 V' 


(16). is then still in agreement with the theoretical value—cf. footnote to (11). 
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theory (and from this analysis an empirical 6’ has no meaning when f 
exceeds 2). If f is small enough that 4/ (1—3/f) can be expanded we 
have for the constant terms of (18) and (19) 


Domb and Salter’s expression for the zero energy, 2R0’, with 6’ the 
constant limiting value, now corresponds to the use of the constant 
limiting @ in (5): 

€y=2RO(1—4f) S 2 RO’. 


Thus if the empirical 226’ is still a good approximation to the zero energy 
it is because the limiting @ is a good approximation to the new Debye 
parameter, and does not indicate that the customary Debye formulae 
retain their validity. 

The new second moment can not be found from a force constant % as 
before since we have no linear wave equations for the w,;; instead, by 
(1), we may write formally 


Lover Sy ae Us : 
[e200 a ON Urs 4s Seb ao: Pern, ph oe (20) 
x clearly depends on the vibrational energy, so that (19) is of the same type 
as the theoretical formula (13) for 0. 

For helium, Domb has replaced U,) in (16) by the energy at 0°K, which 
includes the essential part of the vibrational term y; thus (19) with 
(20) gives the basis for his result in the new theory, provided 6” is replaced 
by the proper 6—and provided we may work from the first few terms 
only of (14) or (17) as a good approximation in this particular case. 
Certainly for solid helium the parameter x of the above analysis is by no 
means small—for the range of 7’ (~ 10°K) to be studied in the following 
paper 0’/T’ ~ 7-20. But Domb and Salter show that for a face-centred 
cubic lattice x may be as large as 6 and yet the first terms of their 
expansion be adequate, a result* depending on relations between the 
higher order moments of the lattice spectrum w,°. We might assume a 
face-centred cubic structure for helium without much thermodynamical 
error, and possibly the moments of the new harmonic spectrum w,; 
(which depend on the fourth degree terms in the lattice potential as well 
as on the terms of second degree) will have similar properties. In that 
case the larger values of ¢;=fw,/k7’ which occur in the helium spectrum 
would be of little effect, even though the exact expression 


k(3é,/sinh 3&5)” 
ee a 
* The radius of convergence of this expansion, or of that quoted below, is 
27; yet no question of convergence really arises since in any case only a few 
terms would be used. Domb and Salter quote for the coefficient of (hw), /k7') 
in (15) a value of the order 10~*, even though I,°, for example, is equal to 4—a 
change in order of magnitude which is very surprising. 
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for the specific heat of the jth oscillator (~ 0 for €; ~ 10 or more) may 
obviously differ for such values from the first few terms of the expansion 


1 1 
ne Salbe eee crates 
/ 1269 + 3q9 6 cn 


which underlies (14). But nothing very useful can be said without 
further knowledge of the new moments p,, and the possibility of Domb 
and Salter’s analysis remains very doubtful in this particular case. 
For helium, then, we must regard (13) as essentially a theoretical formula 
for which the empirical parameter @ (got by fitting the integral expression 
(6) to experimental specific heats) provides an approximate solution. 
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Paramagnetic Resonance in Sodium and Potassium Superoxides 


By J. E. Bennett, D. J. E. Incram and M. C. R. Symons 
University of Southampton 


and P. GrorGcsE and J. STANLEY GRIFFITH 
Department of Colloid Science, Cambridge 


[Received February 3, 1955] 


THE superoxides of sodium and potassium, NaO, and KO,, are of interest 
because the O, ion has an unpaired electron (Neumann 1934 and 
Kassatochkin e¢ al. 1936). Such compounds should give rise to para- 
magnetic resonance absorption, and we have recently detected this in 
both salts. 

The pale yellow sodium peroxide, available commercially, was investi- 
gated first. Chemical and physical tests show that this substance contains 
about 10% of the higher oxide NaO, (George 1942); and it can thus be 
used directly as a magnetically dilute sample of the O, ion, the other 
anions present being diamagnetic. The amorphous mixture was placed 
in a dried cavity resonator, and its absorption studied at 90°k. A broad 
asymmetrical line was obtained, having a maximum intensity at high 
fields, and falling off to a marked shoulder at low fields. Measurements 
were made at 9000 Mc/s, 23 000 Mc/s and 36 000 Me/s, and the separation 
between the low-field shoulder and the high-field maximum was found to 
be accurately proportional to the frequency used. This shows that the 
spread of the absorption line is due to a g-value variation, and not 
an energy level splitting. The g-values obtained were g,=2-175+-0-005 
and g,—2-002+0-005. Within the larger experimental error of --0-05, 
the values obtained from the potassium superoxide, prepared by passing 
oxygen over the heated distilled metal, were identical. Complete 
disappearance of the signal, when the superoxides were treated with 
water, confirmed that it was due to the O, ions. 

We may interpret these results in terms of the electronic structure of 
the ground state of the O, ion. In molecular orbital notation this is 
expected to be 

(c, 1s)? (a, 18)? (oy 28)? (o%, 2s)? (o, 2p)? (7 2p)* (7, 2p)® (Massey 1938). 
As with nitric oxide, the free ion would be in a 7// state, the spin-orbit 
coupling giving a 2/71). and *J/,). state. In nitric oxide *JT,/. lies below 
2TTs).; Og , being one electron short of the closed shell of Fy, will have 
the multiplet inverted with *JJ3. below *J/,j).. However in the solid 
at low temperatures we expect the symmetry about the molecular axis 


2H2 
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to be removed (Zhdanov and Zvonkova 1952), one of the 7, 2p orbitals 
having a higher energy than the other. The magnetic properties of the 
ion will then depend upon the ratio, , of this energy difference to the 
spin-orbit coupling constant. Assuming the ion is not free to rotate, we 
distinguish two limiting cases for the lower state of the ion. For pu zero 
we have the free ion with g,=4, g,=0; for » infinite the orbital contribution 
is completely quenched giving y,=g,—2; and for intermediate values 
2<g<4, 0<g,<2. The detailed theory gives g,=1-99 when g,—2-175, 
in reasonable accord with the results quoted above. Within the ex- 
perimental error, 4. is the same for NaO, as for KO, which is a little 
surprising because they probably differ in crystal structure (Templeton 
and Dauber 1950, Zhdanov and Zvonkova 1952). 

Preliminary measurements on potassium ozonate give a symmetrical 
absorption line with a much smaller g variation, centred on the free-spin 
value. In contrast to this and the superoxides, the absorption obtained 
from frozen solutions of chlorine dioxide is shifted to lower field values, 
and hyperfine components due to interaction with the chlorine nuclei are 
resolved in the wings. This is to be expected as the unpaired electron 
moves over the whole molecule in this case. Further measurements are 
in progress on all these compounds, and the theory and results will be 
reported in detail in the near future. 
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The Decay of 175 and 177 Ytterbium 


‘By H. pr Waarp 
Nobel Institute for Physics, Stockholm* 


[Revised MS. received February 16, 1955] 


THE only spectroscopical measurements on the B--active isotopes of 
ytterbium hitherto reported are those of Cork et al. (1950) on 17°Yb and 
those of McGowan (1951) on 177Yb. These authors have not observed 
the rotational levels found by Temmer and Heydenburg (1954) in lutetium 
by Coulomb excitation. A part of the B-decay from the ground states of 
175 and 177 ytterbium should lead to the first rotational states in the 
lutetium daughter nuclei and y-rays de-exciting these states should 
therefore be found in the decay of these ytterbium isotopes. Such 
y-rays were indeed observed in the present investigation. 

1, 1Yb (4-2d). 99-8% pure Yb,O, was irradiated for one week at 
Harwell. A part of the active material was then further purified by ion 
exchange chromatography. The y-spectrum of this activity was studied 
with a scintillation spectrometer, conversion and £-spectra were taken 
with a double focussing and an intermediate image magnetic spectro- 
meter, respectively. $y and e~~y coincidences and delayed coincidences 
were measured in the intermediate image spectrometer with a y-recording 
Nal(T1)-crystal placed directly behind the source. An oscilloscope+- 
multiplier arrangement resembling that described by Hofstadter and 
McIntyre (1950) was used for recording coincidences between selected 
lines of the y-spectrum and electrons focussed by the magnetic lens on 
the detector Geiger counter. The results obtained from these measure- 
ments are summarized in tables 1 and 2. 


Table 1. y-rays in the Decay of 17°Yb 


y-energy conv. ratios conv. coeff. rel, int. y-rays 
os (kev) 


y, | 113-040-3 | K/(LLM+N)=2540-5 | ap=2:2540-5 | y,/y,=0-540-05 


Yo 281+1 K/L>4 

Ky,/Ky,=0°3 +0:04 Vo/¥v3=0-5 +0-05 
V3 395:140°3 | K/L=5:9+0-7 

L\(M+N)=4-3+1 


The decay scheme of !7°Yb, given in fig. 1(a) could be determined 
unambiguously from the experimental results. The first excited level 
of 175Lu is undoubtedly identical with the first rotational level observed 
by Heydenburg and Temmer and must therefore be assigned 9/2-+-. 
The conversion and intensity data of y, and ys, show that these y-rays 
very probably are both mixtures of #14 M2-radiation (with 


* Communicated by Professor Manne Siegbahn, 
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M2/E1\=442% for y, and M2/H1=26+7% for ys), in which case a 
7/2— or 9/2— assignment must be made for the second level. A measure- 
ment of the angular correlation between y, and ys, described in the 
accompanying letter (Hartmann et al. 1955), show that a 9/2— assignment 
is most probable. 


Table 2. 8-components in the Decay of 17°Yb 


max. en. rel. int. (%) log fT’ in coine. with 
pocoaes (kev) (approx.) 
Bo 463 +3 87 6-25 +0-05 
1 3 7-3 +0-2 V1 
Ba 7243 10 4-8 +0-2 Yo Ys 


M51 y 


(a) (b) 
Disintegration schemes of 175Yb (a) and !77Yb (b). The index numbers of the 
B- and y-rays refer to tables 1 and 2 (17°Yb) and 3 and 4 (177Yb). 


From delayed coincidence measurements performed both with the 
coincidence arrangement mentioned above and with a fast delayed 
coincidence circuit of similar construction as one described earlier 
(de Waard 1952), upper limits of 6 x 10-8 sec for the half life of the 395 kev 
level and of 2 10~® sec for that of the 113 kev level could be determined. 

2. 177Yb (1-8h). This activity was produced by fast neutron irradia- 
tion of the 99-8°%, pure Yb,O, in the big cyclotron of this Institute. 
The activated oxide was used without further purification. The 
y-spectrum of this activity was studied with a scintillation spectrometer, 
conversion and f-spectra were taken in an intermediate image magnetic 
spectrometer. {$y coincidence and delayed coincidence measurements 
were made with scintillation detectors for both B- and y-rays, utilizing 
the oscilloscope coincidence arrangement used also for the measurements 
on '7°Yb, From delayed coincidence plots a half life of 0-122-+-0-005 psec 
was deduced for the state at 146 key in !77Lu, in good agreement with 
the value of 0-13-+0-02 see found by McGowan (1951). Further 
results obtained from the measurements are given in tables 3 and 4, 


ee 
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The decay scheme of !77Yb, given in fig. 1(b), follows from the experi- 
mental results. The first excited state of 177Lu must be the first 
rotational level of that nucleus (9/2+). The conversion data of Vo 
_ Seem to be in accordance only with a 90% H 1+10% M2 character for 
this ray. The 146 kev level should therefore have a negative parity and 
a spin 5/2, 7/2 or 9/2. 


Table 3. y-rays in the Decay of 177Yb 


4 y-energy conv. ratios conv. coeff. rel. int. y-rays 
ee a (kev) (K/L+M) 

V1 119+1 = lj 

Vo piigoal 35h &p—=0-63+0-08 | y,/y2=0-18 +0-04 


Table 4. {$-components in the Decay of 177Yb 


max.en. | rel. int. (%) log fT 


EOIN! (Mev) (approx.) 
B, | 1:30-40-05 88 6-20 40-05 
Be 33 7-55 +0-2 
Bs 9 7-05 £0-2 


From the character of y, and the lifetime of the 146 kev level follows 
a square of nuclear matrix element | M/|?~1 for de-excitation of this 
level by M2 radiation and | M|?~5x10-* for de-excitation by #1 
radiation (| M |? is expressed in ‘ Weisskopf units’). Such very slow 
£1 transitions have also been found by Sunyar (1955). 

A more extensive report on the work summarized here will be published 
in Arkiv for Fysik. 


The author is much indebted to Professor Manne Siegbahn for the 
privilege to work at The Nobel Institute for Physics, to Mr. W. Forsling 
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Angular Correlation of Gamma Rays in *7°Lu 


By L. Axrruinp, B. Hartmann and T. WIEDLING 
Laboratory of Physics, University of Stockholm, Sweden 


[Received February 16, 1955] 


As a complement to a f-spectroscopic investigation of 17°Yb reported 
by De Waard (1955) we have made an investigation of the angular 
correlation of the two cascaded gamma rays of 281 kev and 113 kev 
energies in !75Lu following the B-decay of 17°Yb. 

The investigation was made with two scintillation spectrometers 
coupled in coincidence. One spectrometer accepted only the full energy 
photopeak of the 113 kev gamma ray. The 53 kev x-ray peak of Lu 
was suppressed by shielding the Nal(T1) crystal in the front with 1-3 mm 
tin. The discriminator in the other channel was put just below the 
281 kev photopeak. About 0:7 mm lead was placed in front of this 
Nal(TI) crystal. With the exception of the lines mentioned above and 
the cross-over line at 395 kev no other lines belonging to the decay of 
175Yb were found resolved by the scintillation spectrometers. 

Pile irradiated Yb,O, was chemically purified by ion exchange chroma- 
tography. The correlation measurements were performed on a source 
of Yb(NO 3), dissolved in water. The volume of the solution was about 
6x 10-3 cm’. Coincidences were collected for every fifteen degrees 
between 90° and 180°. 

A least squares fit of the experimental points to the function 
W(0)=1+2A,,Po5, (cos @), corrrected for the finite angular resolution 
of the counters, was made yielding a normalized experimental correlation 
function W(#)=1--(0-202-+0-012)P, (cos @). A small A, coefficient with 
a value of about —0-004 was found, too. 

The ground state of '7°Lu is known to have spin 7/2 (Mack 1950). 
The first excited state of 17°Lu has been identified as a rotational level 
and must therefore be assigned 9/2 (De Waard 1955). Conversion data 
seems to characterize this transition as a mixture of electric quadrupole 
and magnetic dipole radiation, H2/1/1—0-33-+0-10, and further the 
281 kev gamma ray as a mixture M2/H1=(4+2)% so that the 395 kev 
level must probably be assigned 7/2— or 9/2— (De Waard 1955). 

We therefore tried the two following spin sequences: 9/2+9/2—7/2 
and 7/2>9/2>7/2. It is immediately found that pure character of 
both cascaded y-rays is not consistent with the correlation measurement 
(Biedenharn and Rose 1953). Assuming a spin 9/2 for the excited state 
at 395 kev and a pure character of the transition leading from that state 
to the first excited state our correlation measurement gives a mixing 
ratio H2/M1=(17+1)% for the rotational transition to the ground state 
while a spin 7/2 and pure character of the first transition yields 


H2/M1=(3-10-5)%. Both values lie outside the limits of error of the 
conversion measurements, 
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When both transitions are mixed a knowledge of one mixing ratio 
is required to find the mixing ratio of the other (Rose 1954). Using the 
ratio given by De Waard for the second transition the spin sequence 
9/2—+-9/2-+7/2 gives a mixing ratio of 2° for the first transition and 
A,=—0-002. The mixing ratio is here very insensitive to variations in 
the mixing ratio of the 113 kev gamma ray. 

The spin sequence 7/2->9/2->-7/2 gives a mixing ratio of 14% and 
A,=0-012. A comparison with spectroscopic results shows that the 
spin sequence 9/2—+9/2—>7/2 is most probable. 

By making new and more accurate measurements we hope to be able 
to draw more definite conclusions about the spins and transitions of 
175, u. 
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Interaction of Dislocations and F-Centres in Sodium Chloride Single 
Crystals 


By S. Ametinckx, W. Van DER Vorst, R. GEVERS and W. DEKEYSER 


Laboratorium voor Kristalkunde, der Rijksuniversiteit, Rozier, 6, Gent 
(Belgium) 


[Received February 8, 1955] 


Hepass AND MircH ey (1952) have given the first direct evidence for the 
existence of networks of dislocations. They have shown that photolytic 
silver aggregates along certain lines in AgBr single crystals during the 
‘print out ’ process and they identified these lines as dislocation lines. 

We have now succeeded in revealing dislocation lines in the interior 
of melt grown and natural rock salt crystals. This was achieved by 
additive coloration by means of sodium metal and by the precipitation 
of small additions of BaCl,. 

For the coloration of the crystals by sodium metal, we used, in a 
slightly modified form, a technique described by Rexner (1932). 

This method consists essentially in enclosing a certain quantity of 
sodium metal in a cavity made in the crystal. After a suitable heat 
treatment colloidal sodium metal seggregates out along certain surfaces 
and along certain lines within the crystal. 

Photograph 1 (photographs 1-4, Plate 6) is an example of a set of 
very low angle grain boundaries as observed in a well annealed ‘ single 
crystal’ by means of the ultra microscope, at low magnification (100 x ). 
The heavy lines correspond to grain boundaries which are practically 
normal to the plane of observation. The lines of intersection of these inner 
boundaries with the cleavage plane correspond exactly to the pattern of 
etchpits which is revealed after etching with methyl alcohol (Amelinckx 
1954). In many cases the inner surfaces visibly consist of networks of 
lines, as can be seen on photograph 4 (500 x). 

Observations on polygonized specimens prove still more conclusively 
that the observed pattern is due to dislocations. The visible inner 
surfaces (photograph 2, 100 x) have exactly the orientation of the walls 
of dislocations built up during the polygonization of bent specimens. 
Photograph 3 gives at higher magnification (500 x) an example of a 
pattern obtained from such a specimen. The walls have the orientations 
(110) and (110) i.e. normal to active glide planes. Some of them (e.g. AB) 
visibly consist of well separated dots. When focusing lower and lower 
into the crystal, the location of these dots does not change appreciably, 
suggesting that they mark parallel dislocation lines normal to the plane 

- of observation. 

Similar results as described above can be obtained by the addition of a 
small quantity of BaCl, (0-002 to 0-003 g BaCl,.2H,O per 1 g NaCl) 
to the NaCl-melt. On slow cooling, the impurity (BaCl, or a mixed 
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Photograph | 
surfaces in NaCl single crystal revealed by additive coloration with 


sodium metal (100 x). 


Inner 


Photograph 2 
Dislocation walls due to the polygonization of a bent NaCl single erystal (100 x). 
Photograph 3 
Dislocation walls due to polygonization of a bent NaCl single crystal (higher 
magnification, 500 x ). Some of the walls (AB) visibly consist of separated 
lines perpendicular to the plane under observation. 


Photograph 4 


Inner surfaces in deformed and annealed NaCl single crystal (higher magnifica- 


tion, 500 x). 
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LIV. REVIEWS OF BOOKS 


Isotope Geology. By K. Rankama. (London: Pergamon Press Ltd.) [Pp. 535.] 
Price £3 15s. Od. 


‘Isotope GroLocy’ is defined by Rankama as “‘investigation of geological 
phenomena by means of stable and unstable isotopes of elements and of changes 
in their abundance’. Such geological phenomena include determination of 
geological age, generation of heat by radioactivity, measurement of paleo- 
temperatures, origin of coal, calcite etc. The appearance of ‘ Nuclear Geology ° 
(H. Faul, Editor (John Wiley)) close on the heels of ‘ Isotope Geology ’ and 
which covers some of the same field, together with the possibility of a third 
book, all testify to the intense activity which has recently developed in the 
field covered by ‘ Isotope Geology ’. 

Rankama points out that most of the research has in the past been carried 
out by physicists and chemists and the first part of the book (chapters on 
Atoms, Nuclei and Nuclides, Radioactivity, Mass Analysis, Isotope Fraction- 
ation, Properties of Nuclear Radiation) is addressed essentially to geologists. 
The second part (Natural Science of Nuclides) represents the heart of the book. 
Here a chapter is devoted to each element. Literature coverage is unusually 
thorough and in all some 1244 entries appear in the bibliography. One is 
impressed by the enormous disparity in the amount of investigation which has 
been carried out on different elements. Carbon in particular, has been intensely 
investigated and some 330 C1?2/C18 ratios are listed. 

Although addressed primarily fo the geologist, ‘ Isotope Geology ’ will be 
of interest to a much wider audience and all those seeking either general 
information in the field of isotope geology, a thorough reference source, or a 
guide to problems awaiting investigation—and these indeed are many—will 
find the book invaluable. baile he s\.. 


Introduction to Atomic and Nuclear Physics. By HrNry Semat. (London: 
Chapman and Hall Ltd.) Price 50s. net. 


Tuts is a very useful American book for a reader who wants to look up hurriedly 
the evidence for the existence for a neutrino, the properties of the V particle 
or the difference between the various types of particle accelerator. So far as 
the reviewer can see it is reliable and completely up-to-date. Inevitably in a 
book that covers all atomic and nuclear physics in 500 pages the treatment 
is superficial, and at times it is not completely clear; for instance in the 
discussion of fission (p. 404) it is not stated whether the critical energy for fission 
is a quantity that has been observed or calculated. 

As a minor criticism: is it still necessary to include an appendix on the 
evaluation of the integrals [dq which occur in the old quantum theory ? 
And one might wish that an author of a text book would set a good example on 
how to write formulae in a way that makes as little use of hand-setting as pos- 
sible. On page 380 we find 
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which offends against almost every rule that the editors of the Philosophical 
Magazine would commend to their contributors. N. F.M. 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 


